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Abstract 

This paper describes the measurement of the W boson mass, Mw, and decay width, Tw, from the 
direct reconstruction of the invariant mass of its decay products in W pair events collected at a 
mean centre-of-mass energy of y^=172.12 GeV with the OPAL detector at LEP. Measurements of 
the W pair production cross-section, the W decay branching fractions and properties of the W decay 
final states are also described. A total of 120 candidate W^W~ events has been selected for an 
integrated luminosity of 10.36 pb~^. The W"'"W~ production cross-section is measured to be crww = 
12.3±1.3(stat.) ib0.3(syst.) pb, consistent with the Standard Model expectation. The W+W~ — > qq^£ 
and W^W^ qqqq final states are used to obtain a direct measurement of Fw = 1.30lo;55(stat.) ± 
0.18(syst.) GeV. Assuming the Standard Model relation between Mw and Fw, the W boson mass 
is measured to be Mw = 80.32 it 0.30(stat.) it 0.09(syst.) GeV. The event properties of the fully- 
hadronic decays of W"'"W~ events are compared to those of the semi-leptonic decays. At the current 
level of precision there is no evidence for effects of colour reconnection in the observables studied. 
Combining data recorded by OPAL at ^/s ~ 161-172 GeV, the W boson branching fraction to hadrons 
is determined to be 69.8]']3'2(stat.) ib0.7(syst.)%, consistent with the prediction of the Standard Model. 
The combined mass measurement from direct reconstruction and from the W''~W~ production cross- 
sections measured at ~ 161 and ^/s ~ 172 GeV is Mw = 80.35 ± 0.24(stat.) ± 0.07(syst.) GeV. 
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1 Introduction 



In 1996, the LEP collider at CERN entered a new phase of operation, LEP2, with the first e'''e~ 
collisions above the W^W~ production threshold. Approximately 10 pb~^ of integrated luminosity 
was delivered to each of the four LEP experiments at a centre-of-mass energy ^/s ~ 161 GeV. Subse- 
quently, a further ~10 pb^^ was delivered at a higher energy of y/s ~ 172 GeV. This paper describes 
measurements of the production and decay properties of W^W~ events in this higher energy run. 

One of the principal goals of the LEP2 programme is the measurement of the mass of the W 
boson, Mw. Comparison between this direct measurement and the value of Mw determined indirectly 
from precise electroweak analyses at ^/s ~ 91 GeV (LEPl) and elsewhere will eventually provide an 
important new test of the Standard Model (SM) of electroweak interactions. Direct measurements of 
Mw from hadron colliders currently yield 80.33 it 0.15 GeV |jl]-^. Ultimately it is believed that 
LEP2 can achieve a precision on the W mass of approximately 30-40 MeV. The mass dependence 
of the W+W^ production cross-section at threshold, y/s ~ 161 GeV, was used to extract the first 
measurements of Mw from LEP2 data |^-^. At higher centre-of-mass energies measurement of the 
W"'"W~ cross-section itself provides an interesting test of the non-Abelian gauge structure of the 
Standard Model, which predicts substantial destructive interference between the different W'^W" 
production diagrams, thus avoiding unitarity violation at higher energies. 

This paper describes the measurement, above W^W^ production threshold, of Mw and of the W 
decay width, Fw, by direct reconstruction of the invariant mass spectrum of the W decay products. 
The mass and width measurements use events in the channel^ W''~W~ qqqq and W'*"W" qq^f 
(£=e, yu or r). For each event, the mass of the W is reconstructed from its decay products, and Mw 
and Fw are determined by a fit to the resulting distribution. 

The structure of the paper is as follows. Section ^ contains a brief overview of the OPAL detector 
and the Monte Carlo models used. Section ^ describes the selection methods used for the identification 
of samples of W+W" "I/^/, W+W~ qq^£ and W+W" qqqq events at -v/i ~ 172 GeV. 

In Section ^ these selections are used to determine the W^W~ production cross-section and the 
branching fractions of the W boson into various final states. These results may be interpreted in terms 
of a measurement of the Cabibbo-Kobayashi-Maskawa (CKM) mixing parameter |V^s|. In Section |5| 
the W''~W~ — > qqqq and W^W^ qq^f samples are used for a determination of the W mass and 
width. Finally, Section ^ discusses the properties of the hadronic W decays, which are potentially 
sensitive to interesting QCD final-state interactions between the products of the two W decays. 



2 The OPAL Detector, Data and Monte Carlo Models 
2.1 Detector 

A detailed description of the OPAL detector has been presented elsewhere |p and therefore only the 
features relevant to this analysis are summarised here. Charged particle trajectories are reconstructed 
using the cylindrical central tracking detectors, which consist of a silicon microvertex detector, a high 
precision vertex detector, a large volume jet chamber and thin z-chambers. The silicon microvertex 
detector consists of two layers of silicon strip detectors, allowing at least one hit per charged track in 

^Throughout this paper, a reference to W""" or its decay products implicitly includes the charge conjugate states. 
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the angular^ region | cos6\ < 0.93. It is surrounded by a vertex drift chamber. Outside this Hes the 
jet chamber, about 400 cm in length and 185 cm in radius, which provides up to 159 space points 
per track, and measures the ionisation energy loss of charged particles, dE/dx [|^. The z-chambers, 
which improve considerably the measurement of charged tracks in 6, are situated immediately beyond 
and co-axial with the jet chamber. Track finding is nearly 100% efficient within the angular region 
|cos^| < 0.97. The entire central detector is contained within a solenoid which provides an axial 
magnetic field of 0.435 T. 

The electromagnetic calorimeter (ECAL) measures the energy of electrons and photons and pro- 
vides a partial energy measurement for hadrons. It consists of a cylindrical ensemble of 9440 lead 
glass blocks arranged such that the inter-block gaps point slightly away from the origin, and of two 
endcaps, each having 1132 lead glass blocks aligned parallel to the beam axis. The barrel encompasses 
the angular region | cos6\ < 0.82 whilst the endcaps cover the region 0.81 < | cos9\ < 0.98. 

Calorimeters close to the beam axis measure the luminosity using small angle Bhabha scattering 
events and complete the geometrical acceptance down to 34 mrad from the beam axis. These include 
the forward detectors which are lead-scintillator sandwich calorimeters and, at smaller angles, silicon 



tungsten calorimeters [11| located on both sides of the interaction point. 



The iron return yoke of the magnet lies outside the electromagnetic calorimeter, and is instrumented 
with streamer tubes as a hadron calorimeter (HCAL). Muon detectors are situated outside the hadron 
calorimeter. Muons with momenta above 3 GeV usually penetrate to the muon detectors. In addition, 
up to nine hits may be recorded for minimum ionising particles traversing the hadron calorimeter, 
further enhancing muon identification. 



2.2 Data and Monte Carlo 



The basic data selection, luminosity measurement, Monte Carlo (MC) models and detector simulation 
are identical to those described in The accepted integrated luminosity, evaluated using small 
angle Bhabha scattering events observed in the silicon tungsten forward calorimeter, is 10.363 it 
0.045(stat.) ± 0.036 (syst.) pb'^ H, of which approximately 1 pb'^ was collected at 170.3 GeV and 
9.3 pb~^ at 172.3 GeV. The luminosity weighted mean centre-of-mass energy for the data sample is 
^/i = 172.12 ± 0.06 GeV ITl. 



The Gentle |14| semi-analytic program is used to calculate the Standard Model W^W^ cross- 
section which is used throughout this paper to determine the expected number of W^W~ events. The 
use of Gentle is motivated by the fact that it provides a more complete calculation than current 
Monte Carlo generators @. The calculated cross-section is 12.4 pb at ^/s = 172.12 GeV using the 
current world-average W boson mass of Mw = 80.33 GeV ||l],l2|. 



In the analyses described below, a number of Monte Carlo models were used to provide estimates 
of efficiencies and backgrounds as well as the shapes of the W mass distributions. The majority of 
the Monte Carlo samples were generated at ^/s = 171 GeV with Mw = 80.33 GeV. Unless stated 



otherwise, all Monte Carlo samples were generated with a full simulation of the OPAL detector [15|. A 
number of Monte Carlo studies were performed without detector simulation, referred to as generator 
level. 



^The OPAL right-handed coordinate system is defined such that the origin is at the geometric centre of the jet 
chamber, z is parallel to, and has positive sense along, the e~ beam direction, r is the coordinate normal to z, 9 is the 
polar angle with respect to +z and <j) is the azimuthal angle around z. 
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The separation between signal and background processes is complicated by the interference between 
the W'''W~ production diagrams (class^ CC03) and other four-fermion graphs. For example, the 
process TP j"^ — > qq where a is radiated off one of the quarks can interfere with W"'"W~ — > qq^£ and 
W"'"W~ qqqq final states. Monte Carlo samples of W+W" events, restricted to the CC03 diagrams, 
were obtained with the KORALW 



16[, ExCALiBUR |Tj, grc4f Pythia fig] and Herwig 



generators. KORALW was used to determine the efficiencies for W'''W~ events for the selections 
presented in this paper. This sample was generated at \fs = 171 GeV using a value for Mw of 
80.33 GeV. A number of Pythia Monte Carlo samples generated with different values of \fs and Mw 
were used to investigate the sensitivities of the analyses to these parameters. Excalibur was used to 
investigate the sensitivity to Fw- 

The main background process, — > qq, was simulated using Pythia, with Herwig used as 
an alternative to study possible systematic effects. Other backgrounds involving two fermions in the 
final state were studied using KORALZ |21] for e"'"e~ — > n'^fi~, e~^e~ t~^t~ and e'''e~ — > vT/'j'j, 
and Bhwide ||2^ for e'''e~ — > e^e~. Backgrounds from processes with four fermions in the final state 
were evaluated using grc4f, Excalibur and Fermisv p3]. Backgrounds from two-photon processes 



were evaluated using Pythia, Herwig, Phojet Twogen [g5[ and the Vermaseren program 
At least two independent Monte Carlo estimates were available for each category of two-photon and 
four-fermion background. To study the influence of interference effects in the four-fermion final states 
the grc4f and Excalibur Monte Carlo generators were used. In both cases samples were generated 
using the full set of interfering four fermion diagrams. These four-fermion samples were compared to 
samples obtained with the same generator using only the CC03 set of W pair production diagrams. 

Throughout this paper it is assumed that the production cross-section at y/s ~172 GeV for the 
Standard Model Higgs boson, H'^, is negligible. This assumption is valid for a Higgs mass above 
80 GeV. Below this mass, the cross-section becomes significant and the W'*"W~ event selections have 
high efficiencies to select H*^ events, particularly for the process H^Z" qqqq which would result in 
a non-negligible background. 



3 Event Selection 



Event selections have been developed to identify efficiently all Standard Model W^W^ final states 
with low acceptance for background processes. The W"'"W~ event selection consists of three distinct 
parts to select fully leptonic decays, W^W~ l~^iyil'~Vf^i ^ semi-leptonic decays, W~^W~ — > qq^£, 
and fully hadronic decays, W^W^ — > qqqq. To ensure that the selections are mutually exclusive, 
only events failing the W"'"W~ — > i~^i'£i'~T'£' selection are considered as possible W~''W~ — > qq^^ 
candidates and only events failing both the W'''W~ i~^i'£i'~T'£' and W"''W~ qq^£ selections are 
considered as possible W+W^ qqqq candidates. The fully leptonic and semi-leptonic selections 
are separated into the individual lepton types giving ten selected final states: e'^Uee~T'e, e'^VcfJ-'i^fi, 
e+i/gT"!/^, iJi^i/fj_ii~V^, fi~^i>f^T~T'T-, t^Ut-t'Vt-, qqez7c, qq/"i^/x) qq''"'^T and qqqq. No attempt has been 
made to identify the fiavour composition of the hadronic final states. 

The hadronic and semi-leptonic event samples are used to determine the decay width and the mass 
of the W boson and also to study hadronic event properties of W boson decays. All final states are 
used in the measurement of the W pair production cross-sections and decay branching fractions. The 
classification into separate leptonic final states allows the measurement of the W branching ratios with 
or without the assumption of charged current lepton universality. 

^In this paper, the W pair production diagrams, i.e. t-channel Vc exchange and s-channel TP /"f exchange, are referred 
to as "CC03", following the notation of 
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An overview of the event selections is given below. Emphasis is placed on the performance of the 
selections and, in particular, the systematic uncertainties in the selection efficiencies and accepted 
background cross-sections. More detailed descriptions of the W"'"W~ qq^^ and W"'"W~ qqqq 
event selections can be found in the Appendices. 



3.1 Treatment of Systematic Uncertainties 

The KORALW Monte Carlo program (CC03 diagrams only) was used to estimate the efficiencies of 
each of the selections. Two types of systematic uncertainties in the selection efficiencies have been 
considered: generator uncertainties and differences between data and the Monte Carlo simulation, in- 
cluding detector simulation. The generator uncertainties were estimated by comparing the efficiencies 
for four different Monte Carlo generators, KORALW, grc4f, Excalibur and Pythia. The uncertain- 
ties arising from the beam energy dependence of the selection efficiencies were assessed by comparing 
Pythia samples generated at different beam energies. The dominant uncertainties due to the beam 
energy dependence arise because the Monte Carlo samples used to evaluate efficiencies were generated 
at a centre-of-mass energy of 171 GeV rather than 172.12 GeV, the mean energy at which the data 
were recorded. Similarly, the propagation of uncertainties on Mw to the selection efficiencies was 
estimated using samples of Pythia events (including detector simulation) generated with different 
values of Mw- 

For each selection, the four-fermion backgrounds were determined using the grc4f Monte Carlo 
program. The systematic uncertainties on the four-fermion backgrounds were estimated from the 
difference between the accepted cross-sections predicted by the grc4f and Excalibur Monte Carlo 
programs. For both the W"'"W~ — > qqiUi and W"'"W~ — > qqqq selections, the dominant systematic 
uncertainty on the background was due to the modelling of the Z^/j — > qq process. 



3.2 W+W- e+iyd'~17^> Events 



3.2.1 Selection 



Approximately 11% of W pair events are expected to decay through the fully leptonic channel. These 
events may be observed as an acoplanar pair of charged leptons with missing momentum. The selection 
is sensitive to the six possible classes of observed leptons, e+e", /i"*"//", t~^t~, e^//^, e^r^, /i^r^, 
which are expected to be produced in the ratio 1:1:1:2:2:2. The main backgrounds are e''~e~ —>■ Z^Z^, 
e"'"e~ — > Z''e"'"e~, e+e" — > Wel/g, e"'"e~ t^t^ and e"*"e~ e~^e~£^£~ . 

The experimental signature for W^W^ i^uii'^uii events is also the signature for a number 
of non-Standard Model processes. The selection used to identify W^W^ i'^ui^l'^V^i events is 
based on the general charged lepton pair selection used in the search for anomalous production of 
lepton pair events with significant missing transverse momentum |27|. A high efficiency for selecting 



fully leptonic events is obtained by forming the logical "or" of two distinct analyses^. The first 
analysis performs a general selection of events conforming to the acoplanar di-lepton topology. The 
second consists of several selections, each designed to select a particular di-lepton class. This analysis 
was optimised to identify events consistent with being W"'"W~ £'^v^£'~Vii . Both analyses require 
evidence that a pair of charged leptons has been produced in association with an invisible system 

* Referred to as Selection I and Selection II in Reference For Selection I the additional W'''W~ selection criteria 
are applied. 
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that carries away significant transverse momentum. The first analysis simply requires leptonic events 
with significant missing transverse momentum, whereas the second analysis exploits the kinematics 
of W"''W~ — > l^i/ii'^Vii events and the flavour mixture of the expected backgrounds. This leads to a 
complementary acceptance. 

From KORALW Monte Carlo studies, 80% of the expected efficiency for W+W^ £+i/^^'~z7£/ 
events is common to the two analyses, 6% is exclusive to the first analysis, and 14% to the second. 
Both analyses have similar efficiencies for W^W^ i~^i'£i'~'i7£i events with two stable leptons, while 
the second analysis performs significantly better for events with one or two taus. 

3.2.2 W+W- £+ui£'~V£> Classification 

Selected W"'"W~ £'^l'iI'~T'£i candidates are classified as one of the six possible di-lepton combina- 
tions according to the electron and muon identification results, the scaled energy of the leptons and 
the charged track and cluster multiplicity associated with each lepton. Jets are defined using the 
cone algorithm |^]. Those not identified as either electrons or muons are classified as taus. A small 
fraction of selected events are expected to be reconstructed with only one jet. These are events where 
one of the charged leptons is outside the geometric acceptance of the tracking chambers. For these 
one-jet events, the unobserved lepton is taken to be a muon. For two-jet events, corresponding to most 
of the selected events, the classification is further elaborated to improve the assignment of secondary 
electrons and muons from tau decays to the tau class. Identified electrons or muons with energy scaled 
by the beam energy less than 0.3 are classified as taus. In order to recuperate inefficiencies in the 
electron and muon identification algorithm, jets failing the electron and muon identification, but with 
scaled energy greater than 0.5 are reclassified as electrons if E/p exceeds 0.5 and as muons if E/p is less 
than 0.5. Lastly, jets with three or more charged tracks or three or more associated electromagnetic 
clusters (consistent with three-prong tau decays) are classified as taus. 

3.2.3 Results and Systematic Errors 

A total of eight events is selected by the combined event selection at ^/s = 172 GeV as W+W" —>■ 
£~^i'i£'~T'i' candidates. Seven of the events are common to both analyses, while one event is exclusively 
selected by the first analysis. 

Representative kinematic distributions for the selected events are given in Figure [l|, together with 
the Monte Carlo expectations. Figure |l|(a) shows the energy of each charged lepton scaled by the 
beam energy. Figure ||(b) shows cos0_ — cos^+, where 9- and 0+ are the polar angles of negatively 
charged and positively charged leptons respectively. The selected events favour positive values of 
cos0_ — cos^+, as expected for W'''W~ production. These kinematic distributions are consistent with 
the Standard Model expectations. 

The Monte Carlo selection efficiencies for each di-lepton combination are shown in Table [l|. The 
overall efficiency for selecting W^W~ i~^i'ii'~17i' events is determined to be (78.3ib0.4ib2.5)%, where 
the errors are statistical and systematic, respectively. Since events can be rejected on basis of having 
a significant energy deposit in the forward detectors, a correction factor of 0.99 with a systematic 
uncertainty of 0.01 has been applied to the selection efficiencies and the background estimate to 
account for detector occupancy due to ofF-momentum beam particles. Systematic uncertainties on the 
efficiency were estimated based on comparisons of the efficiencies for different Monte Carlo models 
(3.0% for e'^i'ee~T'e, 4.5% for r+z/T-r"!?^ and 2.0% for the other fully leptonic final states). The Monte 



7 



Carlo models have different implementations of both initial and final state radiation effects and the 
modelling of tau decays. In addition, systematic errors were assigned to account for data/Monte Carlo 
agreement (0.8%) and the knowledge of the trigger efficiency (0.4%). 

The expected background cross-sections from Standard Model processes are given in Table ^. The 
systematic errors are based on comparisons of different Monte Carlo generators. Table |^ lists the 
expected and observed numbers of events in each final state. The 8 fully leptonic events observed in 
the data are consistent with the Standard Model expectation. 

3.3 W+W- qqiPi Events 

3.3.1 Selection 

Semi-leptonic final states, W^W^ qq^i^^, are expected to comprise 44% of W^W~ decays. W^W^ - 
qqel^e and W+W" — > qq/xr'^ events are characterised by two well-separated hadronic jets, a high mo- 
mentum lepton and missing momentum due to the unobserved neutrino. The signature for W^W" — > 
qqrT'T- events is two well separated jets from the hadronic W decay and one low multiplicity jet typ- 
ically consisting of one or three tracks. The expected missing momentum is less well defined due to 
the additional neutrino(s) from the decay of the tau. 

The W+W^ — > qqii'e event selection, described in Appendix consists of three separate selec- 
tions, one for each type of semi-leptonic decay. The W'^W" — > qqrl^T- selection is applied only to 
events which fail the W"''W~ — > qqel^e and W^W~ qq/^i^^i selections. Each selection proceeds in 
four stages: 

• Identification of the Candidate Lepton: The track with the highest probability of being a 
lepton from either the decay W — > eUe or W — > fiVfj, is identified. The W"'"W~ — > qqrz7T selection 
uses the track (or tracks) most consistent with being from a tau decay from W tUt-. 

• Preselection: Cuts are applied to the data to reduce the background from qq events. 

• Relative Likeliliood Selection: Relative likelihood selections, based on kinematic variables 
and the lepton candidate, are used to distinguish signal events (W~^W~ — > qqel^e, W"'"W~ 
qq/iI7^ and W^W~ c^TVr) from TP /'^ — > qq background events. Events passing the likelihood 
selection are considered to be W+W~ qq^£ candidates. 

• Event Categorisation: Having suppressed the background using the relative likelihood selec- 
tion, a second relative likelihood is applied to the W"*'W~ — > qqel^e and W+W~ — > qq/^i^/^ can- 
didates in order to categorise them as either W"'"W^ qqel^e, W"''W^ ^qfJ-^fi or W^W~ — ^ 
qqrz/T- events. 

3.3.2 Results and Systematic Errors 

Table |^ shows the efficiencies of the selections for W^W^ — > qq^£ events after categorisation into 
the different channels. These efficiences include corrections which account for observed differences 
between the data and the Monte Carlo simulation. The uncertainties include both systematic and 
statistical contributions. 
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Table |5| lists the sources of the uncertainties evaluated for the selection efficiencies. The respective 
efficiencies of the W"'"W~ — > qqel^e, W+W" c^qfJ-T^^ and W+W" qq7"i^T selections determined 
from Pythia, Koralw, Excalibur and grc4f are found to be consistent within errors. Efficiency 
corrections and systematic errors arising from discrepancies between data and Monte Carlo simulation 
were determined by studying data and Monte Carlo "mixed events" formed by superimposing LEPl 
hadronic decay events and single hemispheres from LEPl events identified as decays to charged 
lepton pairs. Overall corrections to the efficiencies of 0.987ib0.006, 0.999ib0.004 and 1.025ib0.007 were 
found for the W"'"W~ qqel^ej W"'"W~ — > qqfu'^ and W^W~ ^qrur selections respectively. The 
differences between data and Monte Carlo result in a migration of events from the W^W^ qqel^e 
selection to the W'''W~ — > qqTi'T selection. The main difference arises from the simulation of the 
electron identification variables used in the relative likelihoods. The efficiencies summarised in Table 
^ include these corrections. 

Table ^ shows the corrected background cross-sections and total uncertainties for the three selec- 
tions. The systematic errors on the expected background cross-sections are dominated by differences 
between data and Monte Carlo for the two-fermion backgrounds and by differences between generators 
in the case of the four-fermion backgrounds. 

The W^W^ qq^£ selection efficiency for Z^/j — > qq events is small, ~ 0.1%. The background 
estimate from Herwig is consistent with that from Pythia. Since the efficiency is very low, the Monte 
Carlo estimate of the background level from this source is likely to be sensitive to the simulation of 
the tails of distributions. For this reason, the Z^/j — > qq background is estimated from the data. The 
background was found to be 1.2 it 0.5 times that predicted by the Monte Carlo. This was determined 
by fitting the observed likelihood distribution in the region 0.25-0.75, for all data recorded away from 
the Z^ peak (130-140 GeV, 161 GeV and 172 GeV ), with signal and background components where 
the shapes of the respective likelihood distributions are taken from Monte Carlo. The quoted error 
includes a systematic component arising from the variation of the fit region. 

Four-fermion backgrounds were estimated using the grc4f generator. Differences between the 
four-fermion background predictions of the grc4f and Excalibur generators were used to assign 
systematic uncertainties. The W boson pair production cross-section for qqez/g was compared with 
that from the full four-fermion treatment and the difference in the accepted cross-sections taken as 
the effective background. This is important since the background from e'''e~ — > Wez^e can interfere 
with W"'"W~ — > qqeX^g. Generator level (i.e. without detector simulation) studies using Excalibur 
indicated that the contributions from non-CC03 diagrams for the qq/Ul7^ and qqTz/r final states are 
negligible for the experimental acceptance. Two-photon backgrounds are included in the e^e~ff four- 
fermion background. For the two-photon background and each class of four-fermion background at 
least two independent Monte Carlo determinations were used and differences between the predicted 
background cross-sections taken as systematic uncertainties. 

In the 172 GeV data sample 19 W+W" qqeUe events, 16 W"'"W~ qqAi^^/x events and 20 
W"'"W~ — > qqTl7^ events were observed in agreement with the expectation shown in Table ^. Figure |2| 
shows the distribution of lepton energy for accepted events categorised as W'''W~ — > qqel^e and 
W+W^ qq/iI7^. Figure ^ shows four example distributions for selected W"''W~ — > qqTl7^ events. 
The observed distributions for the data are in good agreement with the Monte Carlo expectations. 

3.4 W+W~ — > qqqq Events 



9 



3.4.1 Selection 

Fully hadronic decays, W+W~ qqqq, which are expected to comprise 46% of the total W"'"W~ 
cross-section, are characterised by four energetic, hadronic jets and little missing energy. A selection, 
described in Appendix consisting of preselection cuts and a likelihood analysis is used to separate 
W^W^ qqqq events from the background. The preselection removes events which are likely to be 
from the process — > qq. A relative likelihood for an event being from the process W'''W~ qqqq 
rather than from the dominant background is then estimated using seven kinematic variables. 

For the selection used in the determination of the mass of the W boson and the studies of W'''W~ 
event properties a cut is placed on the value of the relative likelihood. The relative likelihood selection 
is designed to maximise the product of efficiency and purity while limiting possible distortion of the 
W mass spectrum. 

3.4.2 Event Weights 

A modified version of the above relative likelihood is used to give a weight, Wi, to each event reflecting 
the probability that the event originates from a W+W~ — > qqqq decay. These event weights, described 
in more detail in Appendix [B.3| , are used to reduce the statistical uncertainty on the cross-section 
measurement. Given the expected average event weights for both W^W~ qqqq and background 
events, and 1%, the measured W+W~ qqqq cross-section, for an integrated luminosity, L, can 
be expressed as 

T pre ) 

where the summation over i corresponds to summing the weights for the observed events which pass 
the preselection, e^^^ is the efficiency of the preselection for W+W~ — > qqqq and cxj^gj is the accepted 
background cross-section after the preselection. The use of event weights results in a 5% reduction of 
the expected statistical uncertainty on the measured W"^W~ qqqq cross-section. 

3.4.3 Results and Systematic Errors 

The efficiency of the likelihood selection for W^W^ — > qqqq events is estimated from the KORALW 
Monte Carlo simulation to be (79.8 it 0.2 ib 1.2)%, where the errors are statistical and systematic 
respectively. The total expected background cross-section, cibgd) and the contributions from different 
processes are given in Table |^. Also shown is the weighted background cross-section which is the 
expected sum of the event weights for background events passing the preselection. For the event weight 
method, the other parameters necessary to determine the cross-section are the preselection efficiency 
for W"'"W~ —>■ qqqq events, e^g^ = (90.3±0.1)%, and the average event weight for W'^'W" qqqq 
events which pass the preselection, = 0.778 it 0.001, where the errors are statistical. 

The systematic uncertainties evaluated for the W+W^ — > qqqq selection efficiency are given in 
Table |^ This table shows uncertainties on the efficiencies and backgrounds for both the likelihood 
cut selection and for the event weight method. The main systematic uncertainty arises from the QCD 
modelling of both the W"'"W~ — > qqqq process and the dominant background from Z^/'j qq. The 
likelihood function has been re-evaluated using an alternative QCD Monte Carlo model (Herwig) and 
by varying the parameters aq, b, Aqcd and Qo of the Pythia Monte Carlo by it one standard deviation 
about their tuned values . The Monte Carlo distributions used for the likelihood probabilities were 
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rescaled so that the means of the distributions agreed with those observed in the 172 GeV data and 
the analysis repeated. The numbers of LEPl data and Monte Carlo events passing the selection, 
with a modified preselection, were compared. The selection was compared for Monte Carlo events 
generated at different beam energies and with different W masses. Finally the effect of rebinning the 
distributions in the likelihood probability calculation was investigated. The differences observed in 
each case are taken as systematic errors. 

The estimated signal efficiency and expected background cross-section are used to determine the 
expected numbers of signal and background events in Table ^. The observed numbers for the pres- 
election, the likelihood cut and the summed event weight are consistent with the expectations. The 
statistical error on the observed summed event weight is calculated as the square root of the sum of 
the squares of the weights for the events passing the preselection. 



4 W+W Cross-Section and the W Decay Branching Fractions 



4.1 172 GeV Results 



The observed numbers of selected W"*"W~ events have been used to measure the W"*"W~ production 
cross-section and the W decay branching fractions to leptons and hadrons. The measured cross-section 
corresponds to that of W pair production from the CC03 diagrams. Additional diagrams can result 
in the same four-fermion final states as produced in the decays of W''~W~ and can therefore interfere. 
When calculating the expected non-CC03 backgrounds for the various event selections the effects 
of this interference have been included. The expected four-fermion backgrounds, quoted throughout 
this paper, include both contributions from non-CC03 diagrams and the effects of interference. The 
four-fermion backgrounds for each final state are calculated from the difference between the accepted 
four-fermion cross-section including all diagrams, and the accepted CC03 cross-section. For this 
determination the grc4f Monte Carlo was used. The associated systematic uncertainty is estimated 
by comparing the predictions of grc4f and Excalibur. At the current level of statistical precision 
such interference effects are small, but not negligible, for the experimental acceptance. In making the 
measurement. Standard Model expectations for four-fermion processes are used. 

This division into CC03 W"'"W~ events and background events is shown in Table ^, for a W 
mass of 80.33 GeV ||l|, and a centre-of-mass energy of 172.12 GeV. The data are consistent with 
the Monte Carlo expectation. The systematic uncertainties on the expected numbers of signal events 
include contributions from the current errors of ±0.15 GeV |l|,§] on Mw and ±0.03 GeV on the beam 
energy |jl3| (0.9% and 0.2%, respectively). A W+W~ cross-section of 12.4 pb predicted by Gentle 
was used, with a theoretical uncertainty of ±2%. Uncertainties in the selection efficiencies, accepted 



background cross-sections and luminosity have been given in Sections 2.2-3.4. 



The W^W~ cross-section and branching fractions are measured using the information from the 
ten separate channels. For each of the ten final states, i, the probability of obtaining the number of 
observed events is calculated as a function of the W+W~ cross-section, aww; and the W branching 
fractions. A likelihood C is formed from the product of these Poisson probabilities, Vi, of observing 
Ni events for a Monte Carlo prediction of /Xj events: 



For the W''~W — qqqq selection the likelihood is calculated using the summed event weights for the 
99 preselected events with a Gaussian error of ±6.6 corresponding to the square root of the sum of the 
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squared event weights. The expected number of events in each of the ten channels, /Xj, can be expressed 
in terms of the luminosity, the total (CC03) cross-section at ^/s of 172.12 GeV, (Tww(172 GeV), the 
W boson branching fractions, Br(W — > X), the background cross-section in each channel and the 
efficiency matrix for the W^W~ selections, eij. The entries in the matrix, ejj, shown in Table are 
the efficiencies of the event selections i for a W+W^ decay of type j. The off-diagonal elements of the 
matrix determine the acceptances of each selection for other "W^W" decays. 

Three different maximum likelihood fits have been performed. In the first case ctww(172 GeV), 
Br(W — > el^e), Br(W — > ^y^) and Br(W tV-t) are extracted under the assumption that 

Br(W ^ e77e) + Br(W ^ /xI7^) + Br(W ^ tV^) + Br(W ^ qq) = 1, 

giving: 



Br(W - 


eI7e) 


- n 107+0-025 


± 0.004 


Br(W - 




- n n76+°-°2^ 

— U.UrD_QQ]^9 


± 0.003 


Br(W - 


rVr) 


n 1 90+0.032 
— U.iZ0_QQ2g 


± 0.005 



cTww(172 GeV) = 12.5 ± 1.3 ± 0.4 pb, 

where the first uncertainty is statistical and the second systematic. The systematic error includes 
contributions from the uncertainties in the efficiency, background cross-section and luminosity. The 
largest systematic error arises from the uncertainty in the backgrounds in the W+W" selections. The 
correlations between the above measurements are less than 30%. In the second fit the additional 
constraint of charged current lepton universality is imposed^ 

Br(W ^ eI7e) = Br(W ^ nV^,) = Br(W ^ rVr) = Br(W ^ iV^), 

giving 

Br(W^^I7^) = 0.1021°:°^^ ^0.003, 
cjww(172 GeV) = 12.3 ± 1.3 ± 0.4 pb. 

This value for Br(W tu^) implies a value for the hadronic branching fraction, Br(W — > qq), of 

Br(W ^ qq) = 0.694+°;°^^ ± 0.008. 

In the final fit, the branching fractions are fixed to their Standard Model values and the CC03 cross- 
section is determined to be 

(Tww(172 GeV) = 12.3 ± 1.3 ± 0.3 pb. 

The CC03 cross-section at 172.12 GeV depends, albeit weakly, on the W mass. Therefore the measured 
cross-section, (Tww(172 GeV), can be used to obtain a measurement of Mw: 

Mw = 80.5l^:^;°| GeV, 

where the first uncertainty is statistical and the second is systematic. The latter uncertainty includes 
a 2% theoretical component. The result of the fit using the Standard Model branching fractions is 
consistent with that obtained by taking the total number of observed events divided by the luminosity, 
subtracting the total expected background cross-section and dividing by the overall effective selection 

^For the current level of experimental precision, the effect of Br(W tvt) being ~ 0.1% lower Q than Br(W — » eUe) 
and Br(W — > /^F^) has been neglected. 
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efficiency of (77.4 it 1.0)%. This gives (Tww = 12.3lj;3 P^' "^here the error is the combined statistical 
and systematic uncertainty. 

The measured W^W^ production cross-section at ^/s = 172.12 GeV is shown in Figure together 
with the recent OPAL measurement of o"ww at ^/s = 161.3 GeV ||5|. The energy dependence of the 
W^W~ production cross-section, as predicted by the Gentle program for a representative W mass 
of 80.33 GeV [|l|, ||, is also given in the figure. The data are seen to be consistent with the Standard 
Model expectation. 



4.2 Combination with 161 GeV Results 



The branching fraction results from the 172.12 GeV data can be combined with the OPAL data 
recorded at 161.3 GeV |5|. Simultaneous fits to the CC03 cross-section at 161.3 GeV, crww(161 GeV), 
the CC03 cross-section at 172.12 GeV and the W boson branching fractions are performed in the same 
way as presented above for the 172.12 GeV data. The event selections, efficiencies, backgrounds and 
observed numbers of events for the 161.3 GeV data have been described previously 0]. The results of 



this combination are summarised in Table 11 



Within the framework of the Standard Model, the W boson branching fractions depend on the six 
elements of the CKM mixing matrix, |Vjj|, which do not involve the top quark |Q| 

Br(W ^ qq) _ / MMw)^ V- ,,, ,2 
{l-Br(W^qq)) - . 

where as(M-w) is taken to be 0.120 it 0.005. The effect of finite quark masses, m^, on the W branching 
ratios are of order (mq/Mw)^ Q) and are therefore neglected. The result, given in Table 0, for the 
branching fraction Br(W qq) obtained from the fit assuming lepton universality gives: 



E 



iVjjf = 2.22+°J^ _t 0.07. 



j=u,c; j=d,s,b 



This is consistent with a value of 2 which is expected from unitarity. Using the experimental knowledge 

2 2 2 2 2 

of the sum, I Vud I +|Vus| +|Vub| +|Vcd| +|Vcb| = l.OSitO. 01, the above result can be interpreted 

as a measure of iVcsl, the least well determined of these elements: 



iVcsl = 1.08l°:|^±0.03. 



This result is consistent in value with and has a comparable uncertainty to other determinations of 
|Vcs| which do not invoke unitarity ||l|,|30|. 



5 Measurement of the Mass and the Decay Width of the W Boson 



The W boson mass, Mw, and decay width, Fw, are determined from fits to the reconstructed in- 
variant mass spectrum of W pair candidate events. A kinematic fit is employed to improve the mass 
resolution for the W"'"W~ — > qqqq and W^W~^ qq^^ channels. The fully leptonic final state is un- 
derconstrained as it contains at least two neutrinos and is therefore not used. Two different methods 
are used to fit the reconstructed mass distribution. In the first, Mw is determined by performing 
an unbinned likelihood fit which uses a Breit-Wigner function to describe the signal shape. Monte 
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Carlo events are used to correct the fit result for biases introduced by the event selection, detector 
resolution and effects of initial state radiation (ISR). In the second method, a reweighting technique is 
employed to produce Monte Carlo mass spectra, including detector simulation, corresponding to any 
given mass and width. A binned likelihood fit is used to determine Mw and by comparing the 
shape of the reconstructed mass distribution from the data to that from the reweighted Monte Carlo 
spectra. This method takes into account all resolution, acceptance and ISR effects and is free of bias 
as long as these effects are simulated correctly. In addition, the extension from a one-parameter fit 
determining Mw to a two-parameter fit determining both Mw and Fw is straightforward. Therefore, 
the reweighting method (RW) is used to derive the central results of this paper for Mw and Fw- The 
simpler Breit-Wigner fit (BW) provides a valuable cross-check of the RW results and systematics. By 
appropriately choosing the fit normalisations, both methods are made independent of cross-section 
information and are sensitive only to the shape of the reconstructed mass spectrum. 



5.1 Invariant Mass Reconstruction 



For each selected W^W~ event, the masses of the two W bosons can be determined from measured 
invariant masses of the decay products. Incorporating the constraints of energy and momentum 
conservation into a kinematic fit significantly improves the invariant mass resolution. The resolution of 
the kinematic fit is further improved by neglecting the finite W width event-by-event, and constraining 
the masses of the two W boson candidates to be equal, thus yielding a single reconstructed invariant 
mass per event, m^ec- For W^W~ events, m-rec is strongly correlated to the average mass of the two 
W bosons in the event^. Incorporating the measured jet masses into the kinematic fit, rather than 
treating jets as massless, gives an improved mass resolution and yields a better agreement between 
the reconstructed and average masses. Cuts on the fit probability remove poorly reconstructed events 
and reduce background. In addition, there is an ambiguity in the choice of the correct jet-jet pairing 
in W^W~ — > qqqq events which leads to a combinatorial background. 



5.1.1 The W+W — > qqqq Channel 

In each selected W"'"W~ qqqq event, four jets are defined using the Durham algorithm These 
jets are used as input to a five-constraint kinematic fit requiring energy and momentum conservation 
and equality of the two W boson masses. For each event, three kinematic fits are performed, cor- 
responding to the three possible jet-jet pairings, and placed in descending order of the resulting fit 
probabilities Pi, P2, Ps- In Monte Carlo studies, the highest fit probability. Pi, corresponds to the 
correct jet-jet pairing in about 68% of W'^W" — > qqqq events. In approximately 25% of events the 
second highest probability fit, P2, corresponds to the correct combination. The lowest probability fit 
is dominated by incorrect combinations and is not used. Both fit methods require Pi > 0.01. In the 
BW method, the second fit is also included in the reconstructed mass distribution if P2 > 0.01 and if 
P2/P1 > 1/3. This prescription is adopted to optimise the ratio of signal to background. In the RW 
method, for events in which both Pi and P2 exceed 0.01, a two-dimensional distribution, (Mi,M2) 
is used, where Mi is the reconstructed mass corresponding to the combination having fit probability 
Pj. For those events in which only Pi exceeds 0.01, a one-dimensional spectrum is used. The num- 
bers of selected events used in the mass analyses are listed in Table ^ for the fit ranges discussed in 



Sections 5.2 and 5.3 



* At the generator level, the distribution of the average mass of the two W bosons is described by the same Breit-Wigner 
function as the distribution of the two separate W masses, multiplied by a phase space factor. 
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5.1.2 The W+W- 



qqfT'^ Channel 



For the selected W+W"— > qqtPi candidate events = e or //), the reconstructed mass is estimated 
by defining two jets in the hadronic system and then performing a kinematic fit, including the jets and 
lepton, using the same constraints as for the W"'"W~ — > qqqq channel. Since the three-momentum of 
the neutrino is not known, this results in a two-constraint fit. Events with fit probability greater than 
0.001 are retained for the analysis. The probability cut is lower than for the W^W^ — > qqqq channel 
because the semi-leptonic channels have less background. In the case of the W"''W~— > qciTu-r channel, 
the energy of the visible tau decay products is often a poor estimate of the tau lepton energy because 
of the presence of additional neutrinos in the tau decay. To account for this in a kinematic fit, the 
direction of the tau lepton is taken to be that of the identified tau jet and the energy of the tau lepton 
is set to half the beam energy, with an error large enough to cover the kinematically allowed range. 
However, Monte Carlo studies show that the mass resolution of this fit is not significantly better than 
that obtained by dividing the invariant mass of the hadronic system by the ratio of the visible energy 
of the hadronic system to the beam energy. Furthermore, the fit introduces more bias and therefore 
both the BW and the RW method use the reconstructed mass as determined by the scaling method 
for the W'''W~^ qqT^'r channel. To reduce background and remove poorly reconstructed events, 
the BW method requires that the kinematic fit probability be greater than 0.01 for W^W^^ qq'T'i^T 
events. The numbers of selected events in these channels are listed in Table 12. 



5.2 Breit-Wigner Fit to the Reconstructed Invariant Mass Spectrum 

The first method used to determine the W mass involves fitting an analytic function to the distribution 
of the masses obtained from the kinematic fit. Because of biases introduced by event selection, detector 
resolution, ISR, and phase-space effects, the method has been calibrated by performing similar fits to 
Monte Carlo samples generated with known W masses. 



5.2.1 Signal Shape 

A variety of analytic forms to describe the signal and background shapes were investigated. From 
Monte Carlo studies the signal shape is found to be well described up to mrec ~ 84 GeV by a function 
based on a relativistic Breit-Wigner function with fixed width, 

^2 -p2 
("irec-"io)^ + "T'?ecr2' 

where mo and F are characteristic parameters of the signal peak. The fits are limited to the range 
40-84 GeV, where the lower boundary is determined by considerations regarding the background 
normalisation and is discussed below. The width, F, should be regarded as a parameter which embodies 
both the width of the W boson and experimental effects and is fixed to the value predicted by the 
Monte Carlo, 3.8 GeV. 



5.2.2 Background Shape 

Background arises mainly in the W^W~ qqqq channel, both from qq events, and from 

incorrect combinations in W'''W~ qqqq events. Monte Carlo studies show that a quadratic form in 
mrec describes the combinatorial background well up to mj-ec ~ 84 GeV. The qq background 
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has a different shape. It was found that this can be described by the function Ap°'e where 
p = -E'beam ~ ^rec the momentum with which a W pair of mass rrij-ec would have been produced 
and where a and f3 are positive parameters. The shapes and relative normalisation of the background 
contributions are determined from Monte Carlo samples, but their overall normalisation is allowed to 
vary in the fit. Extending the allowed fit range down to 40 GeV improves the determination of the 
background fraction in the peak region, as this low mass region is dominated by background. 



5.2.3 Results of the Breit-Wigner Fit Method 



To extract Mw from the reconstructed mass distribution, an unbinned maximum likelihood fit to 
a relativistic Breit-Wigner plus background function is used. The fit is performed using the recon- 
structed mass distribution from all channels combined. For comparison, it is also performed on the 
W^W~ — > qqqq and W"'"W~^ qq^f channels separately]]. The fit results are given in Table 0, and 
compared with the data in Figure |5|. The expected statistical error has been studied using independent 
subsamples of Monte Carlo events corresponding to the same integrated luminosity as the data. These 
studies demonstrate that the error returned by the fit accurately reflects the r.m.s. spread of the mo 
distribution. 



The correction for bias in the BW method is estimated using Pythia Monte Carlo events corre- 
sponding to various input values of Mw and the beam energy, -Ebeam- Small samples of Monte Carlo 
signal and background events corresponding to the integrated luminosity of the data are processed 
through the same event selection, mass reconstruction, and fitting routines as the data. The mean 
correction is determined by fitting a Gaussian to the resulting distribution of (mg* — m^^) where 
is the W boson mass with which the Monte Carlo sample is generated. This correction is found 
to depend linearly on (-Eboam — i^^'^), as shown in Figure ^. A straight line fit yields: 

W+W- ^ qqqq : m^' - m*^"^^ = -0.363 + 0.0911(Ebeam - 
W+W-^ qqeUe : mg* - m*^"^ = -0.107 + 0.0768(^beam - m^*"') 
Combined : mg* - m^"^ = -0.193 + 0.0779(Sbeam - "i^"") 

with the masses and energies expressed in GeV. These relations are used to correct the fitted mass for 



all biases. The fitted masses, corrections and corrected masses are given in Table 13 



5.3 Fit to the Reconstructed Mass Spectrum using a Reweighting Method 

In the second method, the W boson mass and width are extracted by directly comparing the recon- 
structed mass distribution of the data to mass spectra estimated from fully simulated Monte Carlo 
events corresponding to various values of Mw and Fw- A likelihood fit is used to extract Mw and Fw 
by determining which Monte Carlo spectrum gives the best description of the data. In order to obtain 
the Monte Carlo spectrum for arbitrary values of Mw and Fw, a Monte Carlo reweighting technique 
is employed. 



5.3.1 The Monte Carlo Reweighting Technique 

For a sample of CC03 W pair Monte Carlo events produced with a given mass and width, (M^c, 
T^'^), the detector level mass spectrum corresponding to a different mass and width, (M^™, F^™), 

'^The widths are fixed to different values in the two channels. 
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is estimated by assigning to each event a reweighting factor with which it enters the new spectrum. 
This factor is the ratio of the probabihty that the event would be produced assuming the new values 
of (M-^^jF^™) to the probability that this same event would be produced for the input values 
(M^^,r^^). These production probabilities are given by the W pair production cross-section 

P(Mw, Tw; mi, 7712) oc S>V(Mw, Tw; mi);BW(Mw, Fw; W'2)cro(mi, m2, s), 

where mi and m2 are the generator level masses of the two W bosons produced in the event and (Tq 
is the Born level cross-section for producing a pair of W bosons with masses (772.1,772-2) at a centre-of- 
mass energy ^/s. The values nii and 7722 are distributed according to the running- width, relativistic 
Breit-Wigner function, 

1 #rw 

5>V(Mw,Fw;m) = ^ ^— . 

The reweighting factor for the ith event, fi, is then given by 

The Born level cross-section cancels in the ratio since it depends only on mi, m2 and y/s. 

An estimate of the detector level reconstructed mass spectrum for any (M-^™, F^™) can be obtained 
from one such sample alone with a given input (M^^,F^'^). In general, deviat ions from unity in the 
reweighting factors increase as M-^™ and F^™ differ from iW^^ and F^^*-", respectively. For very 
large fi, the reweighted spectrum becomes sensitive to statistical fluctuations in the input sample. To 
reduce the sensitivity to large reweighting factors, a total of 850 000 events from all W"'"W~ Monte 
Carlo samples with fuh detector simulation are used (Pythia, Herwig, Excalibur, Koralw and 
grc4f Monte Carlo generators). These samples cover a range of input {M^'^ ,Tyf^) values, spanning 
Mw = 78.33 GeV to Mw = 82.33 GeV and widths of 1.7 GeV to 2.5 GeV. The inclusion of ah of 
these samples ensures that the reweighted spectra are reasonably smooth over a large range of (M^™, 
F^™) values. 



5.3.2 The Signal and Background Shapes 

The normalised reweighted mass spectra of the individual Monte Carlo samples are combined into 
a single reweighted spectrum by taking the weighted average, bin-by-bin. As an example of the 
performance of this procedure. Figure compares the reweighted mass distribution for the W"'"W~ 
qqeUi channel for (M^^^, F^"^) = (79.33, 2.047) GeV to the spectrum from the fully-simulated Pythia 
MC sample generated with (M^*^,F^^) = (79.33,2.047) GeV. Good agreement is found. 

The background mass spectra are taken from Monte Carlo and are assumed to be independent of 
Mw and Fw- The background distributions are normalised to the expected number of background 
events. The reweighted signal spectra are then normalised such that the total number of signal plus 
background events corresponds to the observed number of events. This is done separately for each 
channel. For this purpose, the W+W~ qqqq events in which both the first and second highest 
probability fits are used are treated as a separate channel from those W'^W" qqqq events in which 
only the highest probability fit is used. 

Most of the Monte Carlo events have been produced for a beam energy of E^^^ = 85.50 GeV, 
whereas the luminosity weighted average beam energy of the data is E^f^l^ = 86.06 GeV. In order to 
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correct for this, the mass used for reweighting is modified to 

Mnew Tij-new / jTrdata t-iMC \ 

while shifting the resulting spectrum by (-^be*^ — -E'beam)- Since not all the fully simulated samples 
are produced with the same centre-of-mass energy, this is done for each sample separately before 
combining the reweighted spectra. A corresponding shift is made to the background spectra. 



5.3.3 Results of the Reweighting Method 



A binned likelihood is formed from the product of the Poisson probabilities for obtaining the numbers 
of events observed in the data in each bin of mass (and width) using the Monte Carlo expectation 
for (Mw,rw)- The log-likelihood curve is determined separately for each of the different channels. 
Since the channels are statistically independent, the results are combined by adding together these 
likelihood curves. The resulting log-likelihood distribution as a function of the mass and the width 
allows the determination of Mw and T-w and their associated statistical uncertainties. 

To reduce further the sensitivity to statistical fluctuations in the determination of the expected 
number of events, only events with a reconstructed mass in the range 65 GeV < nirec < -E'beam are 
included in the fit. For the two-dimensional mass spectrum of the W"^W~ — > qqqq channel, at least 
one of the two reconstructed masses is required to be greater than 70 GeV and the sum of the two 
masses must exceed 120 GeV. The rather limited mass range is further motivated by the fact that 
events with small reconstructed masses have very little sensitivity to Mw and Fw. 

The results of the fit to Mw and Fw using W^W^ qqqq and W^W^ c^qf^i events are: 

Mw = 80.30 ± 0.27 ± 0.09 GeV, 
Fw = 1.30l°:| ±0.18 GeV, 

where the errors are statistical and systematic, respectively. The correlation coefficient between M-yy 



and Fw is 0.06. The systematic errors are discussed in Section ^.4| . The mass spectra of the individual 
channels together with the fitted reweighted spectra are shown in Figure ^ The likelihood contours 
for the statistical errors are displayed in Figure |^ together with the Standard Model prediction. From 
the shape of the contours and the asymmetric errors on the width, it can be seen that the errors on F-w 
for the present statistics are non-Gaussian. The two standard deviation errors on Fw are l^ ^i GeV. 

A one-parameter fit for the mass is performed by constraining the width using the Standard Model 
relation Q 

Fw = (2.0817 GeV)- ^ 



(80.26 GeV)3' 



giving: 



W+W"^qqqq: Mw = 80.08 ± 0.44 ± 0.14 GeV, 
W+W-^ qqiuc : Mw = 80.53 ± 0.41 ± 0.10 GeV, 
Combined : Mw = 80.32 ± 0.30 ± 0.09 GeV. 



The error on the fitted mass is correlated with the width of the reconstructed mass distribution. 
Consequently, the statistical error on the mass resulting from the one-parameter fit is larger than that 
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resulting from the two-parameter fit. Since this measurement is principahy intended to determine 
Mw as a Standard Model parameter, the result of the one-parameter fit is taken as the central result 
of this paper for the mass of the W boson as determined from the direct reconstruction method. The 
value of Mw obtained from the BW fit is consistent with that obtained using the RW method, thus 
providing a useful cross check of the RW analysis. 



5.4 Systematic Uncertainties 



The main sources of systematic uncertainty on Mw and Fw are studied by changing some feature 
of the analysis. The effect of this can typically be assessed in two ways - either by repeating the fit 
to the reconstructed mass spectrum and observing the change in the fit results, or by determining 
the average of the event-by-event change in the reconstructed mass. The latter method tends to have 
greater statistical sensitivity, but does not readily yield a determination of the error on Fw- For this 
reason, the former method was generally preferred in the RW analysis, and the latter in the BW 
analysis. Despite these differences, the overall systematic errors calculated for the two analyses are 
in good agreement. The following descriptions concentrate on the procedure used to estimate the 



uncertainty in the RW fit method. The estimated errors are summarised in Table 14 



Beam Energy: Uncertainties in the LEP beam energy affect the reconstructed mass spectrum through 
the energy constraints imposed by the kinematic fit. The precision on the beam energy is esti- 
mated to be ±30 MeV [|l^]. The influence of this uncertainty on the fit results has been estimated 
by changing the beam energy in independent samples of signal and background Monte Carlo be- 
fore performing the kinematic fits. In addition, the effect of a possible asymmetry between the 
e"*" and e~ energies is found to be negligible. 

Initial State Radiation: The effect of initial state radiation is to increase systematically the recon- 
structed invariant mass, and consequently the resulting fitted mass, due to biases introduced 
through the energy conservation constraint of the kinematic fit. This bias is taken into account 
in both the RW and BW fit methods to the extent that ISR is accurately modelled in the Monte 
Carlo. The KORALW generator is used to estimate the systematic error associated with the 
incomplete modelling of ISR for the RW fits. Distributions of the mean^ Mw per event are com- 
pared in two samples, one including only first order corrections and one including the full second 
order ISR corrections. From fits to these generator level distributions systematic uncertainties 
of 20 MeV for the mass and 60 MeV for the width are assigned. 

Hadronisation: Sensitivity of the results to the choice of the hadronisation model was studied in two 
ways. In the first approach, a single sample of W^W^ four-fermion final states was hadronised 
separately by both Pythia and Herwig and the fit results compared. By using the same four- 
fermion final states the comparison is sensitive only to differences in hadronisation. Alternatively, 
the mean difference between the reconstructed mass and the generator level average mass is 
determined for the Pythia and Herwig W^W^ samples separately. No significant difference 
is observed in either comparison and the statistical precisions of the tests are taken as the 
systematic errors. Varying the hadronisation model for the background is found to have a 
negligible effect. 

Four-fermion Interference Effects: The Monte Carlo samples used to calibrate the BW method 
and in the reweighting procedure of the RW method include only the CC03 diagrams. In order 

®At the generator level, the average mass of two W bosons is used as an approximation to the results of the five- 
constraint kinematic fit at the detector level. 
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to test the sensitivity of the results to the interference between W^W^ diagrams and other 
four-fcrmion processes, the fit results of a sample generated including the full set of interfering 
four-fermion diagrams are compared to one restricted to the CC03 W"'"W~ diagrams alone. The 
comparison is made using both the grc4f generator and the EXCALIBUR generator. A small mass 
shift was found when using Excalibur, while no such indication was found using grc4f. The 
assigned uncertainties accommodate the results from these two models. 

Detector Effects: The effects of detector mis-calibrations and deficiencies in the Monte Carlo simula- 
tion of the data have been investigated by varying the inputs to the kinematic fit over reasonable 
ranges. The range is determined in each case from a detailed comparison of data and Monte 
Carlo using approximately 1.2 pb~^ of data recorded at -y/s ~ Mz during 1996. The following 
effects are considered: 

• The jet energy scale is determined to be accurate to within 1%. The effect on both data 
and Monte Carlo reconstructed masses is found to be negligible. 

• The calibration of the lepton momentum or energy scale is accurate to within 1%. The 
corresponding changes for the semi-leptonic (e//x) channels are taken as systematic errors. 

• The results of the kinematic fit depend on the covariance matrix of the input parameters, 
which have been determined using Monte Carlo events. Studies of decays indicate 
that the jet energy errors in the data are understood to within 3%. No indication of 
mis- modelling of the errors associated with the jet angles or of the corresponding lepton 
parameters was found. As an estimate of the associated uncertainty, each parameter of the 
covariance matrix was varied within its estimated uncertainty. 

• The algorithm used to correct jet energies and angles combines information from tracks, 
and electromagnetic and hadron calorimeter clusters. Three possible variants in the proce- 
dure are studied. Firstly, the mass reconstruction is repeated with the hadron calorimeter 
removed from the analysis. Secondly, an alternative parametrisation of the jet errors is 
employed. Finally, a different correction algorithm for the combination of track and calori- 
metric information is used. The standard correction procedure combines all tracks and 
clusters into jets and then applies an average correction based on parametrisations of the 
total energy in each detector component contributing to the jet. The alternative algorithm 
first associates tracks with clusters then subtracts the associated charged particle energy 
from the individual clusters prior to jet finding. In all of these checks, the mean and r.m.s. 
of the shifts in the reconstructed masses observed in the data are compatible with the ex- 
pectations from Monte Carlo. The statistical precision of each comparison is taken as the 
associated uncertainty, and the largest of these is included in the total systematic error. 

The errors from all of these checks are added in quadrature to give the total systematic error 
due to uncertainties associated with detector mis-calibrations and mis-modelling. This is the 
dominant experimental systematic uncertainty. 

Fitting Procedure and Background Treatment, RW Fit Method: 

Using 800 fully simulated Monte Carlo subsamples with input values of (Mw,rw) = 
(80.33 GeV, 2.085 GeV), each corresponding to the same integrated luminosity as the data sam- 
ple, it was verified that both the one-parameter and the two-parameter fits have no observed 
bias at the precision studied. In addition, 2000 Monte Carlo test samples, each corresponding 
to the same integrated luminosity as the data, were used to study the reliability of the fit errors, 
the validity of the contour levels, and correlations between the various fitted parameters and 
their associated errors. Test samples were constructed by randomly sampling a parametrisation 
of the reconstructed mass spectrum (including detector simulation) in each channel. 

The possibility of residual bias is further investigated as follows: 
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• To check for any biases in the reweighting procedure that could depend on the input mass 
or width, fits are performed on Pythia Monte Carlo samples with input W masses be- 
tween 78.33 GeV and 82.33 GeV. Only the Pythia samples are used in the reweighting 
procedure, excluding the sample being fitted. An equivalent test is done with the different 
width samples. No significant shifts between the fitted and input values are found in the 
W^W~ — > qqti'i channels and the statistical precision of this test is taken as a systematic 
error. In the W"^W~ qqqq channel there is an indication of a small systematic effect 
which is believed to be an artifact of finite Monte Carlo statistics of the samples employed 
in the reweighting procedure. For values of Mw in the range of interest, 80.33 it 0.50 GeV, 
this bias is ±20 MeV, which is taken as the systematic error. 

• The effect of the finite Monte Carlo statistics in the background samples is estimated by 
comparing the results of repeated fits to a Monte Carlo sample using only fractions of the 
background sample. 

• Variations in the normalisation of the predicted backgrounds, based upon uncertainties in 
the background cross-sections and efficiencies discussed in Section ^, were studied. The 
combined background is varied by 20% in the W^W~ qqqq channel, and 40% in the 
W''"W~ — > c^tUi channels. The observed shifts in the fit results are taken as systematic 
uncertainties. 

• The procedure to correct for the difference between the actual beam energy and the beam 
energy with which the Monte Carlo samples are generated was checked by fitting the value 
for Mw in Pythia sample generated with ^/s = 172 GeV using only the Pythia sam- 
ples generated with -y/s = 171 GeV for the reweighting and background estimation. No 
significant bias is observed. The statistical accuracy of this test is taken as a systematic 
error. 

The errors from all of these checks are added in quadrature to give the total systematic error 
due to uncertainties associated with the fitting procedure and background treatment. 

Fitting Procedure and Background Treatment, BW Fit Method: 

The effects of using different analytic functions to describe the signal and background shapes, of 
fixing various fitted parameters, and of varying the fit range are investigated. In each case, the 
shift in the fit mass observed in the data is compared with the expected shift estimated using 
400 Monte Carlo subsamples, each corresponding to the same integrated luminosity as the data. 
The following variations are considered: 

• A simple non-relativistic Breit-Wigner is used for the signal. 

• The relativistic Breit-Wigner is multiplied by a factor p", where a is a free parameter, and 
p is the centre-of-mass momentum of a pair of particles of mass ?7ircc-0 In this case, the 
data are fitted up to the kinematic limit in m^cc with a fixed to the value determined from 
Monte Carlo. 

• The width parameter, T, is left as a free parameter in the fit. 

• The background is represented by a quadratic in mrec, with parameters determined in the 
fit to the data. 

• The background normalisation is fixed to that expected from Monte Carlo. 

• The relative normalisation of the combinatorial and iP /'^ qq backgrounds is varied by 
±50%. 

• The masses of the combinatorial and — > qq backgrounds are displaced by 1 GeV. 

®The mass distribution would be expected to be modified by a phase-space factor oc p. However, this fails to describe 
the reconstructed Monte Carlo distribution, while a factor p" gives a reasonable fit up to the kinematic limit. 
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In each case, the shift in fitted mass observed in the data is consistent with that estimated from 
the Monte Carlo. The r.m.s. deviation of the differences of the observed shift between the data 
and Monte Carlo is taken as a systematic error. 

Colour Reconnection Effects and Bose-Einstein Correlations: As discussed in Q and refer- 
ences therein, a significant bias to the apparent W mass measured in the W''~W~ — > qqqq 
channel could arise from the effects of colour reconnection and Bose-Einstein correlations be- 
tween the decay products of the W"*" and W~ . Neither of these effects are included in the various 
hadronisation models used. There is no consensus as to the magnitude of such effects, though 
in some models they produce a shift in the measured Mw by as much as -1-100 MeV Q. In 
some models, colour reconnection also affects other properties of W+W" qqqq events, for 
example, the charged particle multiplicity. As will be discussed in Section with the sensitivity 
afforded by the present statistics, no indication of such reconnection effects is observed. In the 
absence of an experimental limit to the biases associated with these phenomena, a systematic 
uncertainty of ±100 MeV is assigned to the W^W~ — > qqqq channel due to these effects. As the 
W^W^ qqqq channel comprises approximately half the total signal sample an uncertainty 
of ±50 MeV is assigned in the combined analysis. One possible test of such effects is to com- 
pare Mw measured in the W"'"W~ — > qqqq channel with that determined from W'''W~— > qq^^. 
Within the present experimental errors the results are compatible. Monte Carlo studies on colour 
reconnection effects indicate that possible effects on the width are typically of the same size as 
on the mass. Thus, a corresponding error of ±50 MeV is assigned to the width as determined 
from a combined fit to all channels. 

5.5 Combination with Cross-section Measurements 

The measurements of Mw from direct reconstruction and from the measurement of the W^W^ pro- 
duction cross-section at y^=172.12 GeV, as presented in this paper, may be combined with the value 
obtained from the W+W^ production cross-section at ^/s=161.3 GeV ||5|, 

Mw = 80.40+°:^^ toil ± 0.03 GeV, 

where the errors are the statistical, systematic and beam energy uncertainties, respectively. The un- 
certainty of the LEP beam energy has been updated since Reference as described in Reference [ |l^ . 
The systematic effects of the three measurements are quite different, and therefore they are combined 
assuming that they are uncorrelated, apart from the uncertainty associated with the LEP beam energy, 
which is taken to be fully correlated. The result obtained is: 

Mw = 80.35 ± 0.24 ± 0.06 ± 0.03 ± 0.03 GeV, 

where the errors are the statistical, systematic, reconnection effects and beam energy uncertainties, 
respectively. 



6 Properties of W Pair Events 

Hadronic data in e'''e~ collisions can be characterised by event shape distributions and inclusive 
observables such as charged particle multiplicities and momentum spectra. In e"'"e~ W"'"W~ events, 
although relatively little data have been collected, it is useful to study such characteristics. In addition 
to tests of Monte Carlo models, measurement of the properties of the hadronic sector of W"*" W~ decays 
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allows the question of colour reconnection to be addressed experimentally. At present there is general 
consensus that observable effects of interactions between the colour singlets during the perturbative 
phase are expected to be small. In contrast, significant interference in the hadronisation process 
is considered to be a real possibility. With the current knowledge of non-perturbative QCD, such 
interference can be estimated only in the context of specific models |1S, 20[ 33-^. 



It has been suggested [34,35] that simple observable quantities, such as the charged multiplicity 
in restricted rapidity intervals, may be sensitive to the effects of colour reconnection. Such effects 
were predicted to be undetectable with data samples corresponding to 10 pb~^. More recently it 
was suggested that the effect on the inclusive charged multiplicity itself may be larger than previously 
considered and that the mean hadronic multiplicity in W^W~ qqqq events, {n'^^), may be as 
much as 10% smaller than twice the hadronic multiplicity in W+W^ — > qq^^ events, {n^^). The 
visible effects of such phenomena are expected to manifest themselves most clearly in low momentum 
regions. Therefore studies of the fragmentation function, i.e. the distribution of the scaled momentum, 

= p/Eheam, are also relevant. 

Most studies of sensitivity to colour reconnection have been estimated within the context of a 
given model, comparing "reconnection" to "no reconnection" scenarios; in general, both the size 
and sign of any changes are strongly dependent upon the model considered. At the expense of a 
reduction in statistical sensitivity, such model dependence can be avoided by comparing directly the 
properties of the hadronic part of W^W^ qq^i^^ events with W^W~ qqqq events. In the 
current study, the inclusive charged multiplicity and the fragmentation function are measured and 
compared for W"'"W~ — > qqqq and W^W~ — > a^tPe events. The quantities A(nch) = {n'^^) — 2{n^^'^) 
and A(xp) = {x^^) — (xp^^*^) are also examined. By way of characterising the global properties of 
W^W^ qqqq events themselves, mean values of the thrust distribution, ((1 — T)^'^), and the 
rapidity distribution relative to the thrust axis, are also measured in this channel. Mean 

values are used for these comparisons due to the relatively small size of the current data sample. 

The models of colour reconnection implementedf^ in the event generators Pythia, Herwig and 
Ariadne ||3^ were used to assess the sensitivities of the quantities above to such effects. Observables 
such as ((1 — T)^'^) and dy^*^]), considered in earlier studies of colour reconnection, had a predicted 
sensitivity that was sufficiently small as to be unobservable at present. The models Ariadne and a 
'colour octet' variant^ of Herwig |^] predict shifts in and {x^) similar in size to the statistical 
uncertainty observed on these quantities in the current data. 



6.1 Correction Procedure 



The distributions of rich, Xp, 1 — T and y are corrected for background contamination using a bin-by- 
bin subtraction of the expected background, based on Monte Carlo estimates. Corrections are then 
applied for finite acceptance and the effects of detector resolution, after which mean values of the 
distributions are calculated. Each observable is evaluated using two samples of Monte Carlo events. 
The first includes full simulation of the OPAL detector and contains only those events which pass 
the cuts applied to the data (detector level). The second does not include initial state radiation or 
detector effects and allows all particles with lifetimes shorter than 3 x lO"^'' s to decay (hadron level). 
Both samples are generated at the same ^/s. Distributions normalised to the number of events at the 
detector and the hadron level are compared to derive bin-by-bin correction factors which are used to 



^"^No retuning of any model was performed in generating events in the various reconnection scenarios. 
Merging of partons to form clusters was performed on a nearest neighbour basis, as a partial emulation of the model 



of Reference 
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correct the observed distributions of Xp, 1 — T and y. 



This correction method is not appropriate for the multiphcity distribution, since resolution and 
acceptance effects cause significant migration of charged tracks between bins. Instead, a matrix 
correction is used to correct for detector resolution effects, followed by a bin-by-bin correction which 
accounts for the effects due to acceptance cuts and residual initial state radiation, as in previous OPAL 
multiplicity studies 



The uncorrected multiplicity and thrust distributions for the W"'"W candidate events before 
background subtraction are illustrated in Figures IC and together with the predictions of Monte 



Carlo events including detector simulation. The background prediction is the sum of all Standard 
Model processes, as described by the models used in Section ^. Good agreement is seen between the 
data and predictions from the models. The simulated W+W^ events do not include colour reconnection 
effects. 



6.2 Systematic Uncertainties 

A number of systematic uncertainties have been considered in the analysis, as summarised in Table 
Each systematic uncertainty is taken as a symmetric error and the total uncertainty is defined by 
adding the individual contributions in quadrature. The dependence of the correction procedure on the 
Monte Carlo model is evaluated by comparing results obtained using Pythia, Koralw or Herwig 
as the W~^W~ signal samples. In each case, the same model is used for the subtraction of the small 
W^W~ background contamination in each channel: W"'"W~ — > qqqq or W^W" — > i'^i'£i'~17ii events 
selected as W"'"W~ — > qqiUi, for example. Uncertainties arising from the selection of charged tracks 
are estimated by varying the track selection cuts and repeating the analysis. The maximum allowed 
values of the distances of closest approach to the interaction region in r-cj) and z are varied from 2 cm 
to 5 cm and from 25 cm to 50 cm, respectively, and the minimum number of hits on tracks is varied 
from 20 to 40. The dependence on charged track quality cuts is the sum in quadrature of these three 
effects. 

To test the dependence of the results on the event selection, alternative selections were considered. 
The W^W^ (^ivi selection was replaced by a cut-based selection similar to that used in ||5|, which 
has a lower efficiency but comparable purity to the likelihood selection. For the W+W~ — > qqqq case, 
the standard likelihood selection is modified by adding a cut on the jet resolution parameter 1/34 such 
that 10% of the selected signal Monte Carlo events are removed from the sample. The entire analysis 
was repeated with the alternative selections. 

The background estimates presented in Section ^ have overall uncertainties of order 25% for the 
qqqq and qqfi^^ channels. These were taken into account by scaling the background by ±25% before 
subtraction from the data. An additional uncertainty was included to allow for the fact that bin 
contents could become negative after background subtraction due to low statistics in the data. 

An uncertainty was evaluated for the Monte Carlo tuning or model dependence of the back- 
ground by shifting the multiplicity distribution by ±1 unit or comparing Pythia with Herwig. These 
gave similar shifts in {rich)- Other background contributions were also varied by considering alterna- 
tive two-photon Monte Carlo samples from Pythia, Herwig or Twogen, by comparing the corrected 
values when treating grc4f W''~W~ or four-fermion events as data, or by neglecting the t'^t~ 
background. The corresponding systematic uncertainty was taken to be the sum in quadrature of the 
uncertainties from; the four-fermion study, the effect of neglecting the t^t^ background and 

the largest shift from the use of any two-photon sample. 
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Since most of the Monte Carlo samples used in the study were generated at ^/s = 171 GeV, the 
analysis was repeated with W''~W~ and Z^/'j samples generated at 172 GeV. The effect of varying 
the signal and background cross-sections over the range expected at these centre-of-mass energies was 
found to be negligible. 

As a further systematic check, the mean values (n^^)' ('^ch^'^)' i^p^^'^)^ i^p^)^ 
were evaluated by applying a correction factor to each of the uncorrected values. This correction is 
the ratio between the Pythia prediction without detector simulation or initial state radiation, to the 
corresponding prediction for the same observable when these two effects are included. The change in 
the corrected value is included as an estimate of the systematic error due to the unfolding process. For 
((1 — T)^^), this is the largest single contribution to the systematic uncertainty, as the event selection 
in the W'''W~ qqqq channel rejects events having two-jet like characteristics which are similar to 
the dominant Z^/'y qq background. Hence, part of the thrust distribution is unmeasured in data. 



6.3 Results 



The mean values of the event properties are as follows, where in each case the first uncertainty is 
statistical and the second systematic. 

(n^^) = 38.3 ±1.1 ±0.6 

{n'^^^") = 18.4 ± 0.9 ± 0.4 

A(nch) = ±1.5 ±2.2 ±0.8 

= (3.22 ±0.13 ±0.08) X 10-2 

{x^^^'') = (3.60 ± 0.20 ± 0.11) X 10-2 

A{xp) = (-0.38 ± 0.23 ± 0.10) x 10"^ 
((l-T)^^) = 0.227 ±0.036 ±0.021 
(|/<i|) = 1.033 ± 0.042 ± 0.025 

The values of ((1 — T)^*^) and agree well with the predictions of models, which are 0.219 and 

1.031 for Pythia, and 0.216 and 1.032 for Herwig, respectively. It is interesting to compare the 
mean charged particle multiplicity in W'''W- events with that in qq events at ^/s = 172 GeV, 



which has a value of approximately 26 |41] 



The difference in mean charged multiplicities in hadronic W decays in qqqq and qq^£ events, 
A(nch), is found to be consistent with zero at the current level of statistical precision. Similarly, 
the measurements of the mean scaled charged particle momenta are consistent in the two chan- 
nels. Figures ^(a) and (b) show the corrected fragmentation functions for the W^W^ qqqq 
and W"^W" qq^£ channels, together with predictions from the Pythia and Herwig models. 
The models are in good agreement within statistical uncertainties in both cases. An alternative mea- 
surement of the mean charged multiplicity may be obtained from the integral of the fragmentation 
function. The values determined in this way are (n^^)=38.2 and {n'^^'") =18.2. Figure |T^(c) shows the 
ratio of the Xp^ distribution to twice the Xp^^^ distribution for low particle momenta, Xp < 0.2. The 
ratio is slightly greater than unity in the high statistics region, in agreement with the positive value 
of A (rich) measured above. 

In summary, the W^W~ event properties presented here are in good agreement with expectations 
of standard QCD models. From studies of reconnection phenomena implemented with the Herwig, 
Ariadne and Pythia models, changes up to approximately one statistical standard deviation in 
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the current data may be expected in and (xp^). Shifts several times larger have also been 

predicted for (n^^) |35]. Defining A(nch) and A(xp) using data alone provides a model independent 
test of possible reconnection effects. The maximum shifts in these variables predicted by the models 
considered are at the level 0.5-1.5 standard deviations for the current data set. There is no indication 
of the effects of colour reconnection on these observables at the current level of statistical precision. 



7 Summary 

W pair events produced in e"'"e~ collisions at 172.12 GeV were analysed. A Monte Carlo based reweight- 
ing method was used to fit the reconstructed mass spectra of the selected hadronic and semi-leptonic 
W pair events to obtain the following values for the mass and decay width of the W boson: 

Mw = 80.30 ± 0.27 ± 0.09 GeV, 
Tw = 1.30l°:| ± 0.18 GeV, 

where the first uncertainty is statistical and the second is systematic. The fitted width is approximately 
one standard deviation below the Standard Model prediction for this value of Mw- Constraining Fw 
to its Standard Model relation gives a measurement of the mass of 

Mw = 80.32 ± 0.30 ± 0.09 GeV. 

These results are consistent with a recent result ||4^. A Breit-Wigner fit to the reconstructed mass 
spectra was used to check the reweighting method for the extraction of Mw, and found to give a 
consistent result. Combining the above value for Mw with the OPAL threshold measurement Q of 
the W boson mass and that derived from the measurement of C7ww(172 GeV), gives 

Mw = 80.35 ± 0.24 ± 0.06 ± 0.03 ± 0.03 GeV, 

where the uncertainties are statistical, systematic, colour reconnection and Bose-Einstein, and beam 
energy, respectively. 

The production cross-section for W^W~ events at this energy was found to be 

cjww(172 GeV) = 12.3 ± 1.3 ± 0.3 pb. 

The decay branching fractions of the W boson were measured assuming only Standard Model W decay 
modes. When combined with the previous OPAL results at -y/s=161.3 GeV the following leptonic 
branching fractions for the W boson are obtained: 



Br(W - 


eI7e) 


- 098+°-°22 

— <J.uyo_QQ2o 


± 0.003, 


Br(W - 




- n73+°-°i9 

— U.UM_ooi7 


± 0.002, 


Br(W - 




- 140+0-030 

— U.i4U„QQ28 


± 0.005. 



If charged current lepton universality is assumed then the leptonic and hadronic branching fractions 
are determined to be 

Br(W^^I7^) = 0.10ll|]:[5iJ ±0.002, 
Br(W^qq) = Omstoml ±0-007, 
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where the errors are 100% anti-correlated. The hadronic branching fraction result gives a direct 
measurement of the sum of squares of kinematically accessible elements of the CKM matrix, 



E 



= 2.22l°J^ ± 0.07. 



i=u,c; j=d,s,b 

This can in turn be interpreted as a measurement of the element |Vcs|, 

iVcsl = 1.08l°:|^ ± 0.03, 

by using measured values of the other |Vij|'^ parameters. The measurement of |Vcs| is consistent in 
value with and has a comparable uncertainty to other determinations which do not invoke unitarity . 
The above cross section and branching fration measurements are consistent with the Standard Model 



expectations and with recent results |30|. 



The predicted effects of colour reconnection and Bose-Einstein phenomena in the fully hadronic 
channel, and also their influence on Mw, are model dependent. A first investigation of these effects 
was performed in data alone by comparing the distribution of the fragmentation function, Xp, and 
mean values of the charged particle multiplicity and Xp for W^W~ qqqq and the non-leptonic 
component of W"'"W~ — > c^qiPg events. The mean values obtained in the two channels are found to 
be consistent, with differences: 

AM = {n'^^)-2{n'^^"') = +1.5 ± 2.2 ± 0.8 



A{xp) = (x^^) - {x^"^^") = (-0.38 ± 0.23 ± 0.10) x lO^^. 



In addition, mean values of rapidity and thrust are determined for W"^W~ qqqq events. All 
measurements are in agreement with predictions of standard QCD models. At the current level of 
statistical precision no evidence for colour reconnection effects was found in the observables studied. 
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Appendices 



A W+W qqiu^ Event Selection 

The selection of each flavour of W^W~ qq^£ events is divided into four distinct stages as discussed 
in Section |3.3| , namely identification of a charged lepton candidate, preselection, relative likelihood 
selection to separate signal from the remaining background, and event categorisation. Details of each of 
these stages are given below. Only events which fail the W'''W~ i~^i'ii'~T'ei selection are considered 
as possible W"'"W~ — > qq^£ candidates. 

A.l Identification of the Candidate Electron and Muon 

In order to maximise efficiency no explicit lepton identification is required. Instead, the track in the 
event which is most consistent with being an electron (muon) from the decay W — > el^e (W — > fjiP^ ) 
is taken to be the candidate lepton. This track is identified using two types of variables: a) lepton 
identification variables, for example, number of hits in the hadron calorimeter or specific energy loss 
in the central tracking chamber; and b) variables representing the probability that the lepton arose 
from a W decay, for example, energy and isolation. These variables are used to calculate, for each 
track, the probability that the track arose from a W — > eug decay, P{e), or from a W — > decay, 
These probabilities are the products of probabilities from the individual variables, which are 
determined using probability density functions obtained from W^W~ Monte Carlo. The track with 
the highest value of P{e) is taken to be the candidate electron in the W^W^ qqel'e selection and 
the track with the highest value of is taken to be the candidate muon in the W^W^ qq/iVfj, 
selection. Each event will yield one candidate electron track and one candidate muon track. 

A. 2 Definition of Variables 

Having selected the most likely electron and muon candidates, variables are constructed which are 
used in the preselection and, subsequently, in the likelihood selection. The following variables are 
used: 

• -Eiepti the energy of the candidate lepton. For electrons this is calculated using the electromag- 
netic calorimeter energy; for muons the track momentum is used, 

• cos^ipmis; the cosiue of the angle between the lepton track and the missing momentum vector, 

• I cos^misl; the cosine of the angle the missing momentum vector makes with the beam axis, 

• -^200; the energy in a cone of 200 mrad around the candidate lepton evaluated using tracks and 
ECAL clusters, 

• P{e) or P{fJ,), the electron or muon identification probability for the candidate lepton track, 

• -Rvis, the visible energy of the event scaled by ^/s, 

• 2/23) where is the value of the jet resolution parameter (Durham scheme at which an 
event is reclassified from i jets to j jets. 
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• SpT) the transverse momentum of the event relative to the beam axis, calculated using tracks, 
ECAL clusters and HCAL clusters, 

• P{s'), the probability from a kinematic fit which estimates the invariant mass of the event, \/s', 

• 9jet, the angle between the lepton candidate and the jet axis of the nearest hadronic jet. 

A. 3 Preselection 

Preselection cuts are applied to remove events which are clearly inconsistent with W^W^ — > qqi'l'^ 
decays. In particular, the preselections are designed to remove most two-photon events and a significant 
fraction of the Z°/7 qq background. Slightly different preselection cuts are applied in the W+W" 
qqeVe and W^W~ — > qq/iI7^ selections. Firstly, events are required to have more than five charged 
tracks and more than seven electromagnetic calorimeter clusters. In addition, the main preselection 
cuts are: -Eicpt> 10.0 GeV, 0.30 < iivis < 1.05, the total energy in the forward luminosity monitors 
< 40 GeV, cos6'ipmis < 0.0 and the energy of the highest energy isolated photon in the event < 
(E'iSR — 10) GeV, where Eisr is the expected energy of an initial state photon for radiative events 
with yfs'^Mio. An isolated photon is defined as an ECAL cluster which is not associated to a track 
and which satisfies the isolation requirement of less than 2.5 GeV of energy in a 200 mrad cone around 
the cluster. Finally, loose cuts are made on the lepton identification probabilities and, in the case of the 
W''~W~ — > qqer'e selection, several cuts are made to reduce the background from converting photons. 
The preselection is approximately 92% efficient for W+W" — > qqei^e and W+W^ — > qq^T'fj, events, 
where about half the loss in efficiency arises from cases where the lepton is outside the experimental 
acceptance for well reconstructed charged tracks. The preselection cuts remove approximately 90% of 
the Z^/j — > qq background. 



A. 4 Relative Likelihood Selection 



For events passing the electron preselection, a relative likelihood method is used to distinguish 
W^W~ qqer'e events from the background of iP/^ — > qq events. After the preselection cuts 
there is still a signal to background ratio of less than 0.1. The likelihoods are based on a set of 
variables, Xi, where the observed values are compared to the expected distributions (obtained from 
Monte Carlo) for W"'"W~ — > qqel/g events from which the probabilities, Pi{xi), are obtained. The 
likelihood, L'^^^'^'^^ is calculated as the product of these probabilities for the individual variables used 
in the analysis. The background likelihood, L^^, is obtained in the same manner using Monte Carlo 
distributions for qq events. The relative likelihood, calculated as: 



^qqc/^c + / X 



where the normalisation factor, /, is the ratio of preselected background to signal cross-sections from 
Monte Carlo. The variables used in the W+W~ qqel^e likelihood are: -Eiept, -£^2005 -P(e), 2/23, Rvis, 
Icos^misl) Y^Pt, cos^ipmis, P{s') and 6jet- The qqfiUfi relative likelihood, C^^'^f^^ is obtained in the 
same manner with a slightly modified set of variables; with P{f^) replacing P{e) and by not using 
6 jet- The inclusion of 0jet in the qqel^e likelihood was observed to improve the signal and background 



separation whereas in the qq/^i^^ selection this was not the case. Figure 13 shows the distributions 
of a number of these variables where each plot shows the combination of the distributions for events 
passing the W^W~ —>■ qqei^e and W+W" — > qqfJ-V^ preselections. 
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Figure |T4|(a) shows the distribution of C"^^'^'^^ for events passing the W"'"W~ — > qqel^e preselection. 
The value of the relative likelihood peaks at around one for W"'"W~ qqeFe events and zero for 
2*^/7 — >■ qq events. Figure ^4|(b) shows the equivalent distribution of C^^i^ for events passing the 
W+W~ qqfiV^ preselection. Events with C"^^'"'^'^ > 0.5 or £1^/^^^ > 0.5 are selected as W+W" 
qqel^e or W+W~ qq/^i^^ candidates, respectively. The combination of preselection and likelihood 
selection rejects about 99.95% of the Z''/7 — > qq background and is approximately 90% efficient for 
W"'"W~ — > qqei^e and W+W^ qq/^^^^i events. 



A. 5 Event Categorisation 

Although the above likelihood selections were optimised to separate W+W~ qqel^e and W^W~ — >■ 
qq//I7^ events from the Z^/j — > qq background, they also select approximately 25% of W~*'W^ qqTX'T 
decays. For this reason events passing the qqer'e selection are re-classified as either W^W^ qqer'e 
or W"'"W~ — > qqrI7,- and events passing the qq^z^/^ selection are re-classified as either qqiJ,u^ or qclrVr- 
This assignment enables separate cross-sections for the three W^W~ qq^i channels to be evaluated 
cleanly. The predominant W^W~ — > qqrl?^ contamination in the qqel^e selection arises from cases 
where the tau lepton decays to an electron or decays into a one prong hadronic final state. Therefore, 
for events identified as W"'"W~ qqei^ei two relative likelihoods, analogous to those described in Sec- 



tion A. 4, are constructed using the same variables as were used in the W"'"W~ — > qqez7e likelihood selec- 
tion. The first relative likelihood attempts to separate W+W" — > qqrI7^ — > qq(eI7e^',-)i^r decays from 
W"^W~ qqel^e decays and the second attempts to separate W+W~ — > qqTl7^ — > qq(7r^n7r°^'r)l^r 
from W^W~ qqez7c- If either of these relative likelihoods is greater than 0.5 the event is re- 
categorised as W"'"W^ qqTT/r- A similar procedure is applied to events passing the W"'"W~ — > qqfu'^ 
selection. The result of the event categorisation is that all events passing the qqel^e and qq/Ul7^ selec- 
tions are categorised as either W+W^ qqel^o W+W^ — > qq^I^^ or W+W^ qqrI7^. 



A. 6 W+W qqri^T- Event Selection 

A relative likelihood selection designed to separate W^W^ qqT^'r from Z^/j qq background is 
applied to events which fail the W+W^ — > qqel^e and W+W^ qq/Ui^/x selections described in Sections 



A.1-A.4. Approximately 25% of all events finally selected as W^W qq^i^r are from the analysis 



given in Sections A.1 - [A.5| , the remaining events being selected as summarised below. The W~''W 



qqTVr event selection proceeds in a similar manner to the W"'"W~ qqeX/g selection described above, 
i.e. 'lepton' identification, preselection and relative likelihood selection. In this case no additional 
event categorisation is performed. The selection has been designed to be sensitive to the four main 
tau decay classes: electron, muon, hadronic one prong and hadronic three prong. Consequently 
four separate selections are applied. The lepton identification of the W+W" qqeX^g selection is 
replaced by the identification of the track most consistent with being from W — > tV^. — > (eI7e^'T-)i^T) 
W — > TUr ifJ'i^ fii^T)''^T and W — > rur — > (Tr^nTr^z^,-)!?,- decays. To be sensitive to three prong tau 
decays, the combination of three tracks which is most consistent with a W — > rI7^ — > {27r^iT^ 1^^)17 r 
decay is also identified. The W"'"W~ qqTl'T selection then proceeds as four preselections, one 
for each of the above cases, and four corresponding likelihood selections. The variables used in the 
likelihood are similar to those used above but include more information about the track (or tracks) 
identified as the tau decay e.g. the invariant mass of all tracks and clusters within a 200 mrad cone 
around the track. Figure |l5| shows a sample of the variables used in the likelihood selections. An event 
with a relative likelihood greater than 0.75 for any one of the four tau likelihoods, is categorised as a 
W^W~ — > qqrr'r event. Events passing more than one of the W^W^ qqTl^r likelihood selections 
enter the final event sample only once. 
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B W+W — ^ qqqq Selection 



The selection of W"'~W~ qqqq events described below is divided into two parts, namely a preselection 
and a likelihood based discriminant formed from the combination of seven observables. A modified 
version of this selection, which uses the same preselection and likelihood observables, is used to obtain 
a weight for each event to be a W"'"W~ — > qqqq event. This variant, described in Section B^, leads 
to a small reduction in the predicted uncertainty this channel contributes to the W^W^ cross-section 
measurement and is therefore used in the analysis of cross-sections and branching fractions. 



B.l Preselection 



Candidate events are required to be classified as hadronic and not be selected by either the 
W"^W~ i~^i'ii'~T'gr or W+W~ — > qqiUg selections. Tracks and calorimeter clusters are combined 
into four jets using the Durham jet-finding algorithm, and the total momentum and energy of 
each of the jets are corrected for double-counting of energy To remove events which are clearly 
inconsistent with W"''W~ — > qqqq events, predominantly radiative Z^/'y events, candidate events are 
required to satisfy the following preselection criteria: 



• V7, the fitted invariant mass of the hadronic system, must be greater than 140 GeV. 

• The energy of the most energetic isolated photon must be less than 0.3^/s. 

• The visible energy of the event must be greater than 0.70y^. 

• 2/34 > 0.003, where is the value of the jet resolution parameter at which an event is reclassified 
from i jets to j jets. 

• Each jet is required to contain at least one charged track. 

• A kinematic fit, which imposes energy and momentum conservation and equality of the W 
masses, is performed on all three possible assignments of jets to W candidates in the event. At 
least one of these combinations must result in a convergent fit. 



The efficiency of these preselection requirements for W^W^ — > qqqq events, e^[^, is 90.3%, whilst 
rejecting 96.6% of the TP qq events. The total background cross-section after the preselection, 
^bad^ is estimated to be 4.5 pb. The preselection accepts 99 events in data. 



B.2 Likelihood Variables 



Events satisyfing the preselection criteria are subjected to a likelihood selection, which discriminates 
between signal and the remaining four-jet-like QCD background. The likelihoods are based on a 
set of seven variables, yi. Probability density distributions are determined for each yi, using both 
simulated W"^W~ qqqq and qq events. Each of these yi distributions are used to construct 

probabilities corresponding to the hypotheses that a given event in data is either W''~W~ qqqq or 
Z*^/7 — > qq. The likelihoods, L''^^^ and L^^, are calculated as the product of these probabilities. The 
relative likelihood discriminant itself, £^5, is defined in terms of these two likelihoods as: 



^qqqq _j_ ^qq ' 
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The value of C'^'^'^'^, between zero and one, is used to discriminate between signal and background 
events. 

The following seven variables, which use the characteristic four-jet-like nature, momentum balance 
and angular structure of W^W~ — > qqqq to distinguish events from the remaining background, are 
used to construct the likelihoods: 



• the logarithm of 1/34, 

• the logarithm of 2/45 , 

• the sphericity of the event, 

• the quantity Jmom, defined as 

J _ jpi +P2-P3- Pi) 

where pi are the jet momenta, and the jets are ordered by energy such that pi is the momentum 
of the most energetic jet, 

• the cosine of the modified Nachtmann-Reiter angle, cosOjy-ji (see [p5[), 

• the cosine of the angle between the two least energetic jets, cos 634, 

• the logarithm of the QCD event weight, qcd420 P^li which is calculated using the tree level matrix 
element for the processes e^e~ — > qqqq, qqgg |^^- This quantity should have large values for 
hadronic Z^/'y decays and smaller values for W"^W~ events. 



The distributions of data passing the preselection criteria for each of these quantities is shown in 
Figure p^, together with the predictions of the Pythia Monte Carlo for W"'"W~ — > qqqq and hadronic 
Z*^/7 decays. Good agreement is seen. The resulting likelihood distribution is given in Figure 17(a). 
A cut value of > 0.2 on the likelihood is used for the measurements of the mass and event properties. 



B.3 W+W — > qqqq Event Weights 



For the W"'"W~ cross-section and W branching fraction measurements the 99 events passing the 
W"*" W~ qqqq preselection are assigned weights, between zero and one, reflecting the probability that 
the event is from W"'"W~ — > qqqq rather than background. The use of event weights results in a 5% 
smaller statistical uncertainty on these measurements^. The weighting method of event classification 
is discussed in |47|. The optimal choice of event weight is the probability that the event arose from 
signal rather than background. In the absence of correlations between the likelihood variables, the 
relative likelihood gives the probability of an event being W~'~W~ — > qqqq. The likelihood variables 



described in Section |B.2| are not uncorrelated. For this reason a linear transformation is applied to 
the likelihood variables, in order to obtain a new set of variables where the off-diagonal elements of 
the covariance matrix are zero. These transformed variables are then used to construct the relative 
likelihood function which is used as an event weight. The transformation is performed using the 
orthogonal matrix which results in transformed variables with diagonal covariance matrices for both 



In the case of the W'''W c^qtue selection the discrimination between signal and background is much greater and 
event weights give a negligible improvement in statistical uncertainty and are therefore not used. 
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the signal and background separately. The event weight, Wi, calculated using these transformed 
variables is defined: 



qqqq 

Wi - 



where iVqqqg and A'qq are the expected numbers of preselected W+W^ qqqq and background 
events respectively. The weights for the preselected events are compared to the Monte Carlo expec- 
tation in Figure [T7|(b). In calculating the event weights, the small background contributions from 
{TP /^){7P /^) qqqq and W'^'W"^ qq^^ are included with the dominant qq background 

in the background covariance matrix and the probability density distributions. 
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w+w- ^ 


Efficiency (%) 




82.4 ± 1.1 ±3.2 
83.8 ±0.8 ±2.2 
76.8 ±0.9 ±2.2 
86.7 ± 1.0 ±2.2 
77.0 ±0.9 ±2.2 
60.6 ± 1.4 ±4.6 



Table 1: Efficiency for selecting each final state, evaluated using the Koralw Monte 

Carlo program, after correcting for detector occupancy. The errors are statistical and systematic 
respectively. 



Source 


Cross-section (fb) 


ff (f = e,^,r,zv,q) 


9 ± 1 


e+e~ff (f = e, /n, r, q) 


32 ± 19 




1.2 ± 0.5 


iiqq [£ = n,T,Ue,iy,^,iyr) 


0.6 ± 0.3 




0.4 ± 0.3 


ei+£i-i^e,Ve, (4 = e,^,r) {h ^ £2) 


18 ± 2 


^l + f^/2~% i^i = ejM;"^) (^1^2 ¥= IJ-t) 


16 ± 13 


Total 


77 ± 24 



Table 2: Expected background cross-sections for different processes, assuming an average centre-of- 
mass energy of 172.12 GeV. The uncertainties include statistical and systematic components. 
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SclecU'd as 


Exp(;cled sigutd 


Expected V)aek. 


Tot al 


C)l")S(U'\'txi 


W+W- ^ 


e+z^ce z^c 


1.3 ±0.1 


0.1 ±0.0 


1.4 ±0.1 


2 


W+W" ^ 




2.6 ±0.1 


0.1 ±0.1 


2.7 ±0.1 


2 


W+W- ^ 




2.3 ±0.1 


0.2 ±0.1 


2.6 ±0.1 


1 


W+W- ^ 




1.4 ±0.1 


0.0 ±0.0 


1.5 ±0.1 





W+W- ^ 




2.2 ± 0.1 


0.1 ± 0.1 


2.3 ± 0.1 


1 


W+W- ^ 




0.9 ±0.1 


0.1 ±0.0 


1.0 ±0.1 


2 


W+W- ^ 


qqel/c 


16.8 ±0.4 


1.2 ±0.4 


18.0 ±0.6 


19 


W+W- ^ 




17.4 ±0.4 


0.6 ±0.1 


17.9 ±0.4 


16 


W+W- ^ 


qqri^r 


13.5 ±0.4 


2.7 ±0.7 


16.2 ±0.8 


20 


W+W- ^ 


qqqq 


41.3 ± 1.5 


13.1 ± 2.0 


54.4 ± 2.5 


54.1 


Combined 


99.7 ±2.6 


18.3 ±2.2 


117.9 ±3.4 


117.1 



Table 3: Observed numbers of candidate events in each W+W- decay channel for an integrated 
luminosity of 10.36 ± 0.06 pb^^ at 172.12 ± 0.06 GeV, together with expected numbers of signal 
and background events, assuming Mw = 80.33 ± 0.15 GeV. The numbers for the W+W- — > qqqq 
channel are the sums of the event weights for the 99 events passing the preselection (see text) . The 
predicted numbers of signal events include systematic uncertainties from the efficiency, himinosity, 
beam energy, W+W^ cross-section and M\y, while the background estimates include selection and 
luminosity uncertainties. The errors on the combined numbers account for correlations. 



Selected as 


W+W qqez^e 


Generated as 
W+W" qqfiv^ 


W+W qqri^r 


W+W- qqe!7c 
W+W- q^llM^fj, 
W+W- ^ qqrVr 


85.1 ±0.9% 
0.2 ±0.1% 
4.7 ±0.5% 


0.1 ± 0.1% 
87.6 ± 0.8% 
5.2 ± 0.5% 


3.9 ± 0.3% 
4.4 ± 0.3% 
61.4 ± 1.2% 



Table 4: Selection efficiencies for the different W+W" —>■ qq^^ channels after event categorisation 
showing the cross contamination in each selection. The efficiencies, based on KoRALW, have been 
corrected for differences between data and Monte Carlo. The errors include both statistical and 
systematic uncertainties. 





Sig 


nal efficiency error 


(%) 




Source of uncertainty 


W+W- ^ qqel/e 


W+W- ^ qq/iI7^ 




W+W 


~ qqrVr 


Statistical 


0.3 


0.3 






0.5 


Comparison of MC models 


0.4 


0.2 






0.2 


Data/Monte Carlo 


0.6 


0.4 






0.7 


Mw dependence (±150 MeV) 


0.1 


0.3 






0.1 


Beam energy dependence 


0.5 


0.5 






0.8 


Total 


0.9 


0.8 


1.2 



Table 5: Sources of uncertainty on the W+W qq£i^e selection efficiencies. 
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Background cross-sections (fb) 




W+W qqez^c 


W+W qq/ii^ft 


W+W ^ qqTi^r 


qqeI7e 


± 25 


1 ± 2 


42 ± 19 


qqff 


1 ± 1 


24 ± 5 


40 ± 4 


e+e-ff 


82 ± 27 


6 ± 2 


31 ± 10 


Z°/7 qq 


29 ± 13 


23 ± 10 


143 ± 61 




6 ± 1 


1 ± 1 


1 ± 1 


Combined 


118 ± 39 


55 ± 12 


257 ± 65 



Table 6: Background cross-sections for the W+W" qq^£ selections in fb. The four-fermion 
backgrounds have been split into final states containing no electron (qqff), where f refers to a fermion 
other than an electron, one electron (qqel^e) and two electrons (e+e^ff). The background from e+e~ff 
final states includes two-photon processes. The Z^/'y qq background includes the uncertainty 
on the background correction factor of 1.2 it 0.5. All errors include both statistical and systematic 
contributions. 





Background crc 
Likelihood cut 

Cbgd 


)ss-section (fb) 
Event weight 

pre 


qq 

qqff(f/e) 
qqe+e" 


1 ± 1 

1237 ± 281 
127 ± 46 
18 ± 3 


1 ± 1 

1142 ± 181 
108 ± 47 
14 ± 3 


Combined 


1380 ± 280 


1260 ± 190 



Table 7: Background cross-sections for the W+W — > qqqq channel in fb, assuming an average 
centre-of-mass energy of 172.12 GeV. The uncertainties include systematic contributions. 









Selection 






Likelihood cut 


Event 


weight 


Source of uncertainty 


Aesig (%) 


Ao-bgd (pb) 


Aie^^^vTs) (%) 


A{al'^^m) (pb) 


Monte Carlo models 


0.30 


0.08 


0.21 


0.07 


QCD parameter variation 


0.93 


0.12 


1.14 


0.12 


Data/Monte Carlo y/s = 172 GeV 


0.70 


0.18 


1.59 


0.02 


Data/Monte Carlo ^/s = 91.2 GeV 




0.11 




0.09 


Beam energy dependence 


0.25 


0.11 


0.34 


0.06 


Mw dependence (±150 MeV) 


0.03 


0.00 


0.17 


0.00 


Binning effects 


0.15 


0.04 


0.32 


0.02 


Total 


1.24 


0.28 


2.03 


0.18 



Table 8: Sources of systematic error on the signal efficiency and expected background cross-section for 
the W+W~ — > qqqq likelihood selection, which is used in the W mass and event properties analyses. 
The systematic errors on the event weight based quantities used for the cross-section measurement 
are also given. The total is the quadrature sum of all uncertainties. 



39 



W+W~ —>■ qqqq 
Selection 


Expected signal 


Expected back. 


Total 


Observed 


Preselection 
Likelihood cut 
Event weight 


53.0 ± 1.5 
46.8 ± 1.2 
41.3 ± 1.5 


46.7 ±2.1 
14.3 ± 2.9 
13.1 ±2.0 


99.7 ±2.6 
61.0 ±3.1 
54.4 ±2.5 


99 
57 

54.1±6.6 



Table 9: Comparison of expected and observed numbers of events in the W+W~ — > qqqq channel, 
based on an integrated luminosity of 10.36 pb~^, Mw = 80.33 GeV and a total W+W~ production 
cross-section of 12.4 pb. Results are given for the preselection cuts, for the likelihood cut selection 
used in the mass and event properties analyses and for the event weight analysis used in the cross- 
section measurement. The errors on the expected numbers of events include contributions arising from 
luminosity, crww and Mw uncertainties. The error on the observed event weight is calculated as the 
square root of the sum of the event weights squared for the 99 events passing the preselection. 



Event 








Efficiencies [%] for W+W" ^ 










selection 


+ 

e v^e Ve 


e^-e/x Z^(u 


+ ^ 

e feT 






r Ut-t Ut 


qqel/e 




qqrz7^ 


qqqq 


+ 

e I'c 


72.3 


0.1 


5.9 


0.0 


0.0 


0.2 


0.0 


0.0 


0.0 


0.0 




1.9 


74.1 


4.6 


0.3 


5.3 


1.0 


0.0 


0.0 


0.0 


0.0 


+ 

e feT Vr 


7.6 


4.3 


62.8 


0.0 


1.3 


9.0 


0.0 


0.0 


0.0 


0.0 




0.0 


1.0 


0.1 


78.3 


6.4 


0.3 


0.0 


0.0 


0.0 


0.0 




0.1 


3.8 


0.4 


7.9 


60.3 


5.9 


0.0 


0.0 


0.0 


0.0 


+ 


0.3 


0.3 


2.9 


0.2 


3.7 


44.3 


0.0 


0.0 


0.0 


0.0 


qqcz7c 


0.0 


0.0 


0.3 


0.0 


0.0 


0.0 


85.1 


0.1 


3.9 


0.0 




0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.2 


87.6 


4.4 


0.1 


qqrI7^ 


0.0 


0.0 


0.0 


0.0 


0.0 


0.3 


4.7 


5.2 


61.4 


0.2 


qqqq 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.4 


70.3 



Table 10: Efficiency matrix, ejj, for the 172 GeV event selections determined using KORALW (CC03) 
Monte Carlo events. Each entry represents the percentage of generated events in decay channel i which 
are accepted by the selection for channel j. The numbers for the qq^^ selections have been corrected 
for differences between data and Monte Carlo. The qqqq numbers are calculated as the preselection 
efficiency multiplied by the average event weight for each event type. The i'^uii'~Vit efficiencies have 
been corrected by a factor 0.99 to account for detector occupancy. 
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Fitted parameter 




Fit assumptions : 
Lepton universality 


SM branching fractions 


Br(W euc) 
Br(W ^ fiu^) 

Br(W ^ TUr) 

Br(W ^ 
Br(W ^ qq) 


0.098to o20 =t 0.003 
0.073t[J:j]}? ± 0.002 
0.140l[J:[]^g ±0.005 

0.690l|]:[]^^ ± 0.007 


0.10ll[];[!J|^ ± 0.002 
0.698j:[];[!^^ ± 0.007 




cjww(161 GeV) [pb] 
aww(172 GeV) [pb] 


3.9io:9 ±0.2 
12.6 ± 1.3 ±0.4 


3.7toi ± 0.2 
12.3 ± 1.3 ±0.4 


3.6toi ± 0.2 
12.3 ± 1.3 ±0.3 



Table 11: Summary of cross-section and branching fraction results from the combined 161.3 GeV 
and 172.12 GeV data. The results from three different fits described in the text are shown. The 
correlations between the branching fraction measurements from the fit without the assumption of 
lepton universality are less than 30%. 



Channel 


Selected events in fit range 
BW method RW method 


W"*"W~ qqqq (1 jet-jet pairing per event) 
W"'"W~ qqqq (2 jet-jet pairings per event) 
W"'"W~^qqeI7o + W+W"^ qq/^i^^i 
W+W~^ qqrVr 


30 22 
12 21 
28 31 
14 17 


Total 


84 91 



Table 12: Numbers of selected events in each channel for use in the W mass determination. For the 
W"'"W~ qqqq process, up to two jet-jet pairings are used in each event. For the BW method, only 
events with a reconstructed mass of 40 GeV < nirec < 84 GeV are included, while for the RW method, 
reconstructed masses from 65 GeV up to the kinematic limit are used. 





W+W qqqq 


W+W"^ qqiPi 


combined 


Fitted mass (GeV) 
MC correction (GeV) 
Corrected mass (GeV) 


80.19 ±0.39 
-0.19 ±0.01 
80.01 ±0.43 ±0.14 


80.80 ± 0.37 
-0.32 ±0.01 
80.48 ±0.40 ±0.12 


80.50 ±0.27 
-0.26 ±0.01 
80.24 ±0.30 ±0.10 



Table 13: Summary of fit results and Monte Carlo corrections to Mw for the BW fit method. The 
errors on the fitted mass and corrections are statistical only, while for the corrected mass the statistical 



and systematic errors are both listed (see Section |5.4| for details) 
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BW-fit to Mw 


RW-fit to Mw 


RW-fit to Mw and Fw 


Systematic errors 


fixed width 


SM constrained width 






(MeV) 


qqqq 




comb. 


qqqq 




comb. 


Mw 


Fw 


Beam Energy 


27 


27 


27 


30 


30 


30 


30 


15 


ISR 


22 


22 


22 


20 


20 


20 


20 


61 


Hadronisation 


16 


16 


12 


52 


26 


20 


21 


52 


Four-fermion 


42 


37 


37 


42 


48 


42 


35 


70 


Detector Effects 


67 


87 


54 


62 


76 


48 


48 


98 


Fit procedure 


















and Background 


33 


65 


36 


32 


28 


25 


25 


85 


Colour Reconnection 


















and Bose-Einstein 


100 





50 


100 





50 


50 


50 


Total systematic error 


137 


121 


96 


144 


104 


94 


92 


175 



Table 14: Summary of the systematic errors for the various fit methods. For the one-parameter fits 

to determine Mw, uncertainties are given for the fits to W'''W~ qqqq, W'''W~ qq^£ and for the 
combined sample. For the two-parameter fit, the uncertainties for Mw and Fw are listed separately 
only for the fits to the combined sample. 



Systematic variation 




("ch ) 


A (rich) 


(4^) 
xlO^ 


xlO^ 


A{xp) 
xlO^ 


((1 - r)^^) 




W'*'W~ model dependence 


0.30 


0.16 


0.19 


0.018 


0.035 


0.027 


0.005 


0.006 


Track quality cuts 


0.28 


0.18 


0.26 


0.043 


0.055 


0.065 




0.011 


Event selection 


0.15 


0.25 


0.52 


0.047 


0.067 


0.020 


0.008 


0.017 


Background scaling 


0.05 


0.06 


0.15 


0.018 


0.005 


0.054 


0.008 


0.011 


Z/7 background 


0.27 


0.03 


0.22 


0.012 


0.005 


0.017 


0.007 


0.005 


Other backgrounds 


0.06 


0.07 


0.11 


0.017 


0.009 


0.022 


0.001 


0.002 


Beam energy dependence 


0.15 


0.01 


0.18 


0.001 


0.015 


0.016 


0.004 


0.002 


Unfolding procedure 


0.15 


0.22 


0.29 


0.025 


0.051 


0.026 


0.015 


0.001 


Total 


0.56 


0.42 


0.75 


0.076 


0.11 


0.10 


0.021 


0.025 



Table 15: Individual contributions to the systematic uncertainties on the average event properties. 
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Figure 1: Distributions of kinematic variables for selected W^W~ i^i^gi' V^' events, showing: 
(a) the lepton energy divided by the beam energy, and (b) cos 6- — cos 6^ . The data are shown as 
points with error bars. The Monte Carlo prediction for the sum of W+W" and all other Standard 
Model processes is shown as the open histogram, while the non-W"'"W~ processes arc represented by a 
doubly hatched histogram. The Monte Carlo samples have been normalised to the collected integrated 
luminosity of 10.36 pb~^. 



43 



18 



S 16 



S 14 

CM 
© 

© 12 



10 



8 



2 




n — I — I — I — I — I — T" 



n — I — I — I — I — T" 



OPAL 



10 20 30 40 50 60 70 80 90 100 

Lepton energy (GeV) 



Figure 2: The energy of the lepton for events selected as W"''W~ — ^ qqer'c or W'^'W" — >■ qq/zi^/i 
after categorisation. For electrons this is determined using the ECAL energy and for muons using the 
track momentum. The open histogram shows the Monte Carlo prediction (Pythia) for an integrated 
luminosity of 10.36 pb~^ and total W+W~ cross-section of 12.4 pb. The contribution from W+W~ 
qqrl^T- decays is shown as the single hatched histogram and the contribution from background processes 
as the doubly hatched histogram. The data are shown as the points with error bars. 
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OPAL 
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Lpj (GeV) Hadronic energy (GeV) 

Figure 3: Distributions of kinematic variables for selected W"'"W~ qqTi^T events, showing: (a) 
estimated visible energy of the tau decay, excluding the energy of the neutrino(s), (b) visible energy 
of the event divided by ^/s R^is, (c) net transverse momentum J2pt of the event, and (d) uncorrected 
energy of the two hadronic jets in the event, calculated using tracks and unassociated ECAL clusters, 
when the jet containing the tau candidate is excluded (this variable is not used in the likelihood). The 
open histogram shows the Monte Carlo prediction (Pythia) for an integrated luminosity of 10.36 pb~^ 
and total W+W" cross-section of 12.4 pb. The contribution from W+W" qqeX^g and W''"W~ 
qq/il/ju decays is shown by the single hatched histogram and the contribution from background sources 
as the doubly hatched histogram. The data are shown as the points with error bars. 
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Figure 4: The dependence of (Tww on -/i, as predicted by Gentle for Mw = 80.33 GeV. The W+W~ 
cross-sections measured at ^fs = 172.12 GeV (this paper) and at ^/s = 161.3 GeV are shown. The 
error bars include statistical and systematic contributions. 
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Figure 5: Reconstructed mass distributions for the data. The main plot shows the W^W^ qqqq and 
W"'"W~— > qqiPt samples combined, and the insets present them separately. The solid curves display 
the results of unbinned maximum likelihood fits to a relativistic Breit-Wigner signal plus background 
in the range 40-84 GeV as described in the text. The background function alone is shown by the dark 
shaded region. 
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OPAL Monte Carlo 




Figure 6: Mean of the difference between fitted and true mass from fits to many Monte Carlo sub- 
samples, for Pythia Monte Carlo events generated with various input W masses and at different 
beam energies, plotted as a function of the difference between beam energy and true mass. The 
W^W^ — > qqqq and W+W^^ qq-^z^£ samples are shown separately, as well as both combined. The 
lines represent fits to the points, with parameters given in the text. The arrows indicate the actual 
corrections applied. 
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Figure 7: Illustration of the reweighting procedure for the reconstructed mass spectrum of the 
W"'"W~ — > qqfz7^ signal events. The spectrum of a sample produced with M\y = 79.33 GeV 
is compared with the reweighted spectrum obtained from samples with W boson masses between 
Mw = 78.33 GeV and Mw = 82.33 GeV. Only samples produced with the same Monte Carlo gener- 
ator (Pythia) are used for this comparison and the Myy = 79.33 GeV sample has not been included 
in the reweighting procedure. The Monte Carlo statistics for the final fit including all samples, will 
be larger by more than a factor of 3. 
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Figure 8: Fits to the mass spectra of the different channels using the reweighting method. The dots 
represent the data, the Ughtly shaded histograms the result from the fit for (Mw,rw) and the darkly 
shaded ones the background contributions to it. For the W^W^ — > qqqq candidates in which two 
jet-jet pairings per event are used, the underlying histogram represents the fit result. For clarity, the 
bin size used for the one-dimensional histograms is twice that used for the fit. 
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Figure 9: The 39% and 86% contour levels of the two-parameter fit using the reweighting method. 
The projections of these contours onto the axes give the one and two standard deviation errors quoted 
in the text. The straight dashed line gives the dependence of the width on the mass according to the 
Standard Model. The one-parameter fit for the mass alone is constrained to this line. The solid line 
is the Standard Model prediction for variations of (Mt,op,-MHiggs), both given in GeV. For small fitted 
values of Fw, the 86% contour is very sensitive to the particular distribution of masses observed in 
the data. 



51 



|12 

1 10 
o 

I 8 



1 1 — I 1 — I 1 — I 1 1 1 1 1 1 1 1 1 — I 1 — I 1 — I — I 1 1 1 1 1 1 r 



(a) 



OPAL 







>16 



^14 
d 



8 
6 
4 
2 





J I I L 



10 



60 

for qqqq 




+ OPAL data 

PYTHIA 

KORALW 
HERWIG 
Background 



Figure 10: Uncorrected charged multiplicity distributions for (a) W+W~ — > qqqq events and (b) the 
hadronic part of W+W~ qq^£ events. The points indicate the data, the histograms show the total 
expected signal and background contribution for various signal models, and the hatched histogram 
shows the expected background. The bins are two units wide. 
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Figure 11: Uncorrected thrust distribution for W''"W~ qqqq events. The points indicate the data, 
the histograms show the total expected signal and background contribution for various signal models, 
and the hatched histogram shows the expected background. 
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Figure 12: Corrected Xp distributions for (a) W+W~ qqqq events, (b) the hadronic part 
of W^W" — qq^^ events and (c) the ratio of the W+W~ — > qqqq distribution to twice the 
W"'"W~ — > qq^£ distribution. The points indicate the data, with statistical and systematic un- 
certainties added in quadrature, and the predictions of various Monte Carlo models (without colour 
reconnection) are shown as histograms. c- < 
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Figure 13: Distributions of some of the variables used in the W+W~ — > qqel^e and W+W" — > qq//!^^ 
selections for events passing either the W''"W~ — qqel^e or the W+W" — > qq^V^ preselection cuts. 

The contribution from W^W~ qqTi^T decays is shown as the single hatched histogram and the 
contribution from background processes as the doubly hatched histogram. The data are shown as the 
points with error bars. 
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Figure 14: Values of the W+W~ qqeX/g and W+W" — > qq/Ul^^ relative likelihoods, C°i^^^^ and 
£mt^v^^ for events passing the W+W" — > qqez7e and W+W" — qq/xl^/i preselection cuts respectively. 
The contribution from W^W" qqTi'r decays is shown as the single hatched histogram and the con- 
tribution from background processes as the doubly hatched histogram. The peaks in the background 
distributions near one arise predominantly from four-fermion background. The data are shown as 
the points with error bars. Events having relative likelihood values greater than 0.5 are selected, as 
indicated shown by the arrows. 
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Figure 15: Distributions of some of the variables used in the W"'"W~ qqTX^r selections for events 
passing the W"'"W~ — > qqTz7^ preselection cuts. Some events appear more than once in the plot since 
they can satisfy more than one of the W+W" — ^ qqTi^r preselections. The variables shown are; (a) 
the energy of the tau decay products, excluding the neutrino(s), (b) the energy within a 200 mrad 
cone about the tau candidate track(s) calculated using tracks only, (c) the scaled visible energy, after 
preselection cuts of 0.3 < i?vis < 1-0 and (d) the transverse momentum of the event. The contribution 
from W+W" — qqel/e and W+W" — qqiiV^ decays is shown as the single hatched histogram and 
the contribution from background processes as the doubly hatched histogram. The data are shown as 
the points with error bars. 
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Figure 16: Distributions of the seven variables that have been used in the W"''W~ qqqq hkehhood 
selection, after preselection. The points indicate the data, the open histogram shows the W+W" — 
qqqq Monte Carlo and the hatched histogram shows the total background expectation. 
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Figure 17: Distributions of (a) the likelihood distribution, used to select events for the mass and 
event properties analyses and (b) the event weights, used for the cross-section measurement. The 
points indicate the data, the open histogram shows the W''~W~ — ^ qqqq Monte Carlo and the hatched 
histogram shows the total background expectation. 
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